Jammed Hard-Particle Packings: From Kepler to Bernal 
and Beyond 

S. TorquatcQ 

Department of Chemistry, Department of Physics, Princeton Center for Theoretical 
Science, Princeton Institute for the Science and Technology of Materials, 
and Program in Applied and Computational Mathematics, Princeton University, 
Princeton New Jersey, 08544 USA and 

School of Natural Sciences, Institute of Advanced Study, Princeton New Jersey, 
08540 USA 

F. H. Stillingefl 

Department of Chemistry, Princeton University, Princeton New Jersey, 
08544 USA 



1 



Abstract 

Understanding the characteristics of jammed particle packings provides basic in- 
sights into the structure and bulk properties of crystals, glasses, and granular media, 
and into selected aspects of biological systems. This review describes the diversity 
of jammed configurations attainable by frictionless convex nonoverlapping (hard) 
particles in Euclidean spaces and for that purpose it stresses individual-packing 
geometric analysis. A fundamental feature of that diversity is the necessity to 
classify individual jammed configurations according to whether they are locally, 
collectively, or strictly jammed. Each of these categories contains a multitude of 
jammed configurations spanning a wide and (in the large system limit) continuous 
range of intensive properties, including packing fraction 0, mean contact number 
Z, and several scalar order metrics. Application of these analytical tools to spheres 
in three dimensions (an analog to the venerable Ising model) covers a myriad of 
jammed states, including maximally dense packings (as Kepler conjectured), low- 
density strictly-jammed tunneled crystals, and a substantial family of amorphous 
packings. With respect to the last of these, the current approach displaces the 
historically prominent but ambiguous idea of "random close packing" (RCP) with 
the precise concept of "maximally random jamming" (MRJ). Both laboratory pro- 
cedures and numerical simulation protocols can, and frequently have been, used 
for creation of ensembles of jammed states. But while the resulting distributions 
of intensive properties may individually approach narrow distributions in the large 
system limit, the distinguishing varieties of possible operational details in these 
procedures and protocols lead to substantial variability among the resulting distri- 
butions, some examples of which arc presented here. This review also covers recent 
advances in understanding jammed packings of polydisperse sphere mixtures, as 
well as convex nonspherical particles, e.g., ellipsoids, "superballs" , and polyhedra. 
Because of their relevance to error-correcting codes and information theory, sphere 
packings in high-dimensional Euclidean spaces have been included as well. We also 
make some remarks about packings in (curved) non-Euclidean spaces. In closing 
this review, several basic open questions for future research to consider have been 
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identified. 



*Elcctronic address: torquato@princcton.edu 
^Electronic address: fhs@princeton.edu 



Contents 

I. Introduction 
II. Preliminaries and Definitions 

III. Lessons From Disordered Jammed Packings of Spheres 

IV. Jamming Categories, Isostaticity, and Polytopes 

A. Jamming Categories 

B. Isostaticity 

C. Polytope Picture of Configuration Space 

V. Order Metrics 

A. Specific Order Metrics 

B. Characteristics of a Good Order Metric 

VI. Order Maps and Optimal Packings 

A. Strict Jamming 

B. Collective and Local Jamming 

C. Broader Applications to Other Condensed States of Matter 

VII. Protocol Bias, Loss of Ergodicity, and Nonuniqueness of Jammed States 
VIII. Attributes of the Maximally Random Jammed State 
IX. Packings of Spheres With a Size Distribution 

X. Packings of Nonspherical Particles 

A. Ellipsoid Packings 

B. Supcrball Packings 

C. Polyhedron Packings 

D. Additional Remarks 

XL Packing Spheres in High-Dimensional Euclidean Spaces 
XII. Remarks on Packing Problems in Non-Euclidean Spaces 



a 



4 



XIII. Challenges and Open Questions 



82 



Acknowledgments 



References 



84 



84 



I. INTRODUCTION 



The importance of packing hard particles into various kinds of vessels and the questions it 
raises have an ancient history. Bernal has remarked that "heaps (close-packed arrangements 
of particles) were the first things that were ever measure d in the form of basketfuls of grain 
for the purpose of trading or the collection of taxes" (IBernall . Il965l ). Although packing 
problems are easy to pose, they are notoriously difficult to solve rigorously. In 1611, Kepler 
was asked: What is the densest way to stack equal-sized cannon balls? His solution, known 
as "Kepler's conjecture, " was the fac e-centered-cubic (fee) arrangement (the way your green 



grocer stacks oranges) 



Gauss 



1831 



proved that this is the densest Bravais lattice packing 
(defined below). But almost four centuries passed before Hales proved the general conjecture 
that there is no other arrangement of spher es in three-d imen sional Euclidean space w hose 
density can exceed that of the fee packing ( Halesl . 2005 ); see 



Aste and Weaire 



2008 



for a 



popular account of the proof. Even the proof of the densest packing of congruent (identical) 
circl es in the plane, the two-dimensiona l analo g of Kepler's problem, appeared only 70 years 



ago (jConway and Sloane 



1998 



Roger si . Il964l ); see Fig. [TJ 



Packing problems are ubiquitous and arise in a variety of applications. These exist in the 
transportation, packaging, agricultural and communication industries. Furthermore, they 
have been studied to help understand the symmetry, structure an d macroscopic physical 



prop er ties of c onden s ed ma t ter phases, liquids, glass e s and c rystals (Ashcroft and Mermin . 



1976 



Bernall. 



1960. 



Maver and Maver 



Weeks et al 



1971 



194ol: 



1965 



Chaikin and Lubenskvl. 



Speedy 



1991 



1995 



Sti 



Woodcock and Angelll. 



analysis of heterogeneous materials (ITorquat 



Russel et al. 



1989 



Torquatol . 120091 ). and granular media ( lEdwardsl . Il994j ) . Understanding 



linger et al. 



198 
itcl . 



1961 



Hansen and McDonald!. 



Stillinger and Salsburg . 



1986 



1969 



). P acking pro blems are also relevant 



2002|), colloids] 



Chaikin and Lubenskv 



or the 



1995 



the symmetries and other mathematical characteristics of the densest sphere packings in vari- 
ous spaces and dimensions is a challenging area of long-standing interest in discrete geometry 
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FIG. 1 A portion of the densest packing of identical circles in the plane, with centers lying at the sites of a 
triangular lattice. The fraction of M. 2 covered by the interior of the circles is <f> = n/\/l2 = 0.906899 . . .. The 
first claim of a proof was made by Thuc in 1892. However, it is generally b elie ved that the first complete 



error-free proof was produced only in 1940 by Fejcs Toth; see 
for the history of this problem. 



Rogers 



1964 and 



Conway and Sloane . 



1998 



and number theory ( 



as coding theory ( 



Cohn and Elkies 



Cohn and Kumar 



2003 



2007b 



Conway and Sloane 



1998 



Conway and Sloane 



1998 



. logers 



Shannon 



1964) as well 



1948|). 



It is appropriate to mention that packing issues also arise in numerous biological contexts, 
spanning a wide s pectrum of length scales. This includes "crowdi ng" of macromolecules 



within living ce 



2008 



Torquato 



Is 



Ellis 



lotaxis) ( Nisoli et al 



20011 ). the packing of cells to form tissue (IGevertz and Torquato 



20021 ) , the fascinating spiral patterns seen in plant shoots and flowers (phyl 



2010 



Prusinkiewicz and Lindenmayer 



19901 ) and the competitive set- 



tlement of ter ritories by animals, the patterns of which can be modeled as random sequen- 



tial packings (ITanemura and Hasegawa 



1980 



Torquato! . 120021 ). Figure [2] pictorially depicts 



macromolecular crowding and a familiar phyllotactic pattern. 

We will call a packing a large collection of nonoverlapping (i.e., hard) particles in either a 
finite-sized container or in <i-dimensional Euclidean space M. d . The packing fraction is the 
fraction of space covered by (interior to) the hard particles. "Jammed" packings are those 
particle configurations in which each particle is in contact with its nearest neighbors in such 
a way that mechanical stability of a specific type is conferred to the packing (see Section HVl) . 
Jammed packings and their properties have received considerable attention in the literature, 
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FIG. 2 Biological examples of packing geometries. Left panel: Schematic representation of the approximate 
numbers, shapes and density of packing of macromolecules inside a cell of Escherichia coli. (Illustration by 
David S. Goodscll, the Scripps Research Institute.) Molecular crowding occurs because of excluded volume 
effects due to the mutual impenetrability of all of the macromolecules. This nonspecific steric repulsion is 
always present, regardless of any other attractive or repulsive interactions tha t might occ ur between the 



macromolecules. Thus, molecular crowding is essentiall y a packing problem ( Ellis 



Close-up of a daisy capitulum (flower head), as given by 



Prusinkiewicz and Lindenmaver 



200lh. R ight panel: 



1990. The most 



prominent feature is two sets of spirals, one turning clockwise and the o ther counterclockwise. There are 



sever al models that relate such phyllotactic patterns to packings problems (jPrusinkiewicz and Lindenmaver 



1990) 



both experimental and theoretical. Within the domains of analytical theory and computer 



approach ( 


Bernal . 


1960 




1965 


lid wards. 


1994 


Edwards and Grinev 


. 2 


001; 


Gao et al. 


2006; 


Liu and Naeel 




1998 




Makse and Kurchan 




2002 




O'Hern et al. 




2003 


Parisi and Zamponi. 


2010 




Silbert et al 




2002 


• 


200^ 




Song et al. 


■ 


200S 




Wvart et al. 


200, 


5), which for a given 



packing procedure aims to understand typical configurations and their fre quency of oc 



currence. The other, more recent 



20044 



2001 



2007a. 



2007 



2004d: 



Kansal et al. 



Torquato et al 



y, is th e "geometric-structure' 



2002b 



Torquato et al. 



2003 



approach ( 



Donev et al 



Torquato and Stillinger 



20001 ). which emphasizes quantitative characterization of single- 
packing configurations, without regard to their occurrence frequency in the algorithmic 
method used to produce them. Our primary objective is to review the latter approach, 
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while at the same time to show that these two approaches are complementary in that they 
represent different aspects of the larger context of hard-particle jamming phenomena. In 
particular, a wide range of jammed packing ensembles can be created by the choice of the 
generating algorithm, and the geometric-structure approach analyzes and classifies individ- 
ual members of those ensembles, whether they be crystalline or amorphous at any achievable 
packing fraction 0. 

The process of cooling an initially hot liquid ultimately to absolute zero temperature 
provides a close and useful analogy for the subject of hard-particle jamming. Figure |3] sum- 
marizes this analogy in schematic form, showing typical paths for different isobaric cooling 
rates in the temperature-volume plane. While these paths are essentially reproducible for a 
given cooling schedule, i.e., giving a narrow distribution of results, that distribution depends 
sensitively on the specific cooling schedule, or protocol, that has been used. A very rapid 
quench that starts with a hot liquid well above its freezing temperature will avoid crystal 
nucleation, producing finally a glassy solid at absolute zero temperature. A somewhat slower 
quench from the same initial condition can also avoid nucleation, but will yield at its T = 
endpoint a glassy solid with lower volume and potential energy. An infinitesimal cooling 
rate in principle will follow a thermodynamically reversible path of equilibrium states, will 
permit nucleation, and will display a volume discontinuity due to the first order freezing 
transition on its way to attaining the structurally perfect crystal ground state. By analogy, 
for hard-particle systems compression qualitatively plays the same role as decreasing the 
temperature in an atomic or molecular system. Thus it is well known that compressing 
a monodisperse hard-sphere fluid very slowly leads to a first-order freezing transition, and 
the resulting crystal p hase corresponds to the closest packing arrangement of those spheres 
(IMau and Husd . Il999l ). The resulting hard sphere stacking variants are configurational im- 
ages of mechanically stable structures exhibited for example by the venerable Lennard- Jones 
model system. By contrast, rapid com pression rates applied to a hard -sphere fluid will create 
random amorphous jammed packings (iRintoul and Torquatol . Il996bl ). the densities of which 
can be controlled by the compression rate utilized (see Fig. [TUin Sec. IV below). In both 
the cases of cooling liquid glass-formers, and of compressing monodisperse hard spheres, it 
is valuable to be able to analyze the individual many-particle configurations that emerge 
from the respective protocols. 

We begin this review, after introducing relevant terminology, by specifically considering 
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glass 
transition 
temperature 
range 



Temperature 



FIG. 3 Three isobaric (constant pressure) cooling paths by which a typical liquid may solidify, represented 
in a volume vs. temperature diagram. An infinitesimal cooling rate from the high-temperature liquid traces 
out the thermodynamic equilibrium path (shown in green), including a discontinuity resulting from the 
first-order freezing transition. This reversible path leads to the ground state defect-free crystalline structure 
in the T — > limit. Very slow but finite cooling rate (not shown) can involve crystal nucleation but typically 
creates defective crystals. More rapid cooling of the- liquid (blue curves) can avoid crystal nucleation, passing 
instead through a glass transition temperature range, and resulting in metastablc glassy solids at absolute 
zero. The volumes, energies, and other characteristics of those glasses vary with the specific cooling rate 
employed in their production. 



packings of frictionless, identical spheres in the absence of gravity, which represents an ide- 
alization of the laboratory situation for investigations of jammed packings; see Section [TTTl 
This simplification follows that tradition in condensed matter science to exploit idealized 
mo dels, such as the Ising model, which is regarded as one of the pillars of statistical mechan- 



ics (jDomb 



1960 



Gallavotti . 



1999 



Onsagerl . Il944l ). In that tradition, this idealization offers 



the opportunity to obtain fundamental as well as practical insights, and to uncover unifying 
concepts that describe a broad range of phenomena. The stripped-down hard-sphere "Ising 
model" for jammed packings (i.e., jammed, frictionless, identical spheres in the absence of 
gravity) embodies the primary attributes of real packings while simultaneously generating 
fascinating mathematical challenges. The geometric-structure approach to analyzing indi- 
vidual jammed states produced by this model covers not only the maximally dense packings 
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(e.g., Kepler's conjecture) and amorphous "Bernal" packings, but an unbounded collection 
of other jammed configurations. This approach naturally leads to the inevitable conclusion 
that there is great diversity in the types of attainable jammed packings with varying degrees 
of order, mechanical stability, and density. 

Important insights arise when jammed sphere packings are placed in a broader context 
that includes jammed states of noncongruent spheres as well as nonspherical objects, as 
discussed in Sees. [IX] and |X] These extensions include polydisperse spheres, ellipsoids, 
"superballs," and polyhedra in three dimensions. In addition, this broader context also 



involves sphere packings in Euclidean space with high dimensions (Sec . IXip. which is relevan t 



to error correcting codes and information theory (jConway and Sloand . 119981 : iShannonl . 119481 ). 
and packings in non-Euclidean spaces (Sec. IXIip . Finally, in Sec. IXHH we identify a number 
of basic open questions for future research. 



II. PRELIMINARIES AND DEFINITIONS 



Some basic definitions concerning packings are given here. A packing P is a collection of 
nonoverlapping solid objects or particles in d- dimensional Euclidean space M d . Packings can 
be defined in other spaces (e.g., hyperbolic spaces and compact spaces, such as the surface 
of a d-dimensional sphere), but our primary focus in this review is M. d . A saturated packing 
is one in which there is no space available to add another particle of the same kind to the 
packing. 



FIG. 4 Symmetries associated with three particle shapes. Chords pass through each particle 
centroid. Left panel: A "superdisk" is centrally symmetric and possesses two equivalent principal 
axes. Middle panel: An ellipse is centrally symmetric but does not possess two equivalent principal 
axes. Right panel: A triangle is not centrally symmetric. 
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We will see subsequently that whether a particle possesses central symmetry plays a fun- 
damental role in determining its dense packing characteristics. A d- dimensional particle is 
centrally symmetric if it has a center C that bisects every chord through C connecting any 
two boundary points of the particle, i.e., the center is a point of inversion symmetry. Ex- 
amples of centrally symmetric particles in M. d are spheres, ellipsoids and superballs (defined 
in Section |X]). A triangle and tetrahedron are examples of non-centrally symmetric two- 
and three-dimensional particles, respectively. Figure H] depicts examples of centrally and 
non-centrally symmetric two-dimensional particles. A d-dimensional centrally symmetric 
particle for d > 2 is said to possess d equivalent principal (orthogonal) axes (directions) 
associated with the moment of inertia tensor if those directions are two-fold rotational sym- 
metry axes such that the d chords along those directions and connecting the respective pair 
of particle-boundary points are equal. (For d = 2, the two-fold (out-of-plane) rotation along 
an orthogonal axis brings the shape to itself, implying the rotation axis is a "mirror im- 
age" axis.) Whereas a d- dimensional superball has d equivalent directions, a <i-dimensional 
ellipsoid generally does not (see Fig. HJ). 

A lattice A in M. d is a subgroup consisting of the integer linear combinations of vectors 
that constitute a basis for M, d . In the physical sciences and engineering, this is referred to 
as a Bravais lattice. Unless otherwise stated, the term "lattice" will refer here to a Bravais 
lattice only. A lattice packing Pl is one in which the centroids of the nonoverlapping identical 
particles are located at the points of A, and all particles have a common orientation. The 
set of lattice packings is a subset of all possible packings in M d . In a lattice packing, the 
space M. d can be geometrically divided into identical regions F called fundamental cells, each 
of which contains the centroid of just one particle. Thus, the density of a lattice packing is 
given by 

(i) 



r Vol(F) ' 

where V\ is the volume of a single d- dimensional particle and Vol(F) is the rf-dimensional 
volume of the fundamental cell. For example, the volume vi(R) of a d- dimensional spherical 
particle of radius R is given explicitly by 

, . 7l d ' 2 R d 

where V(x) is the Euler gamma function. Figure [5] depicts lattice packings of congruent 
spheres and congruent nonspherical particles. 
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FIG. 5 Examples of lattice packings (i.e., Bravais lattices) depicted in two dimensions. Left 
panel: A portion of a lattice packing of congruent spheres. Each fundamental cell (depicted as a 
rhombus here) has exactly one assigned sphere center. Right panel: A portion of a lattice packing 
of congruent nonspherical particles. Each fundamental cell has exactly one particle centroid. Each 
particle in the packing must have the same orientation. 

A more general notion than a lattice packing is a periodic packing. A periodic packing 
of congruent particles is obtained by placing a fixed configuration of N particles (where 
N > 1) with arbitrary nonoverlapping orientations in one fundamental cell of a lattice A, 
which is then periodically replicated without overlaps. Thus, the packing is still periodic 
under translations by A, but the N particles can occur anywhere in the chosen fundamental 
cell subject to the overall nonoverlap condition. The packing density of a periodic packing 
is given by 

Nvi 

pvu (3) 



r Vol(F) 

where p = N/Vo\(F) is the number density, i.e., the number of particles per unit volume. 
Figure [6] depicts a periodic non-lattice packing of congruent spheres and congruent nonspher- 
ical particles. Note that the particle orientations within a fundamental cell in the latter case 
are generally not identical to one another. 

Consider any discrete set of points with position vectors X = {r±, T2, . . .} in M. d . Associ- 
ated with each point G X is its Voronoi cell, Vor(rj), which is denned to be the region of 
space no farther from the point at than to any other point rj in the set, i.e., 

Vor(rj) = {r : |r — < |r — i\,|for all Tj G X}. (4) 

The Voronoi cells are convex polyhedra whose interiors are disjoint, but share common faces, 
and therefore the union of all of the polyhedra is the whole of M. d . This partition of space 
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FIG. 6 Examples of periodic non-lattice packings depicted in two dimensions. Left panel: A 
portion of a periodic non-lattice packing of congruent spheres. The fundamental cell contains 
multiple spheres located anywhere within the cell subject to the nonoverlap constraint. Right panel: 
A portion of a periodic non-lattice packing of congruent nonspherical particles. The fundamental 
cell contains multiple nonspherical particles with arbitrary positions and orientations within the 
cell subject to the nonoverlap constraint. 

is called the Voronoi tessellation. While the Voronoi polyhedra of a lattice are congruent 
(identical) to one another, the Voronoi polyhedra of a non-Bravais lattice are not identical to 
one another. Attached to each vertex of a Voronoi polyhedron is a Delaunay cell, which can 
be denned as the convex hull of the Voronoi-cell centroids nearest to it, and these Delaunay 
cells also tile space. Very often, the Delaunay tessella tion is a trianqu lation of space, i.e., 
it is a partitioning of W 1 into <i-dimensional simplices (jTorquatol . 120021 ) . Geometrically the 
Voronoi and Delaunay tessellations are dual to each other. The contact network is only 
defined for a packing in which a subset of the particles form interparticle contacts. For 
example, when the set of points X defines the centers of spheres in a sphere packing, the 
network of interparticle contacts forms the contact network of the packing by associating 
with every sphere a "node" for each contact point and edges that connect all of the nodes. As 
we will see in Sec. IIVI the contact network is crucial to determining the rigidity properties of 



the packing and corresponds to a subclass of the class of fascinatin 



frameworks, namely strut frameworks; see 



Connelly and Whiteley 



objec ts called tensegrity 



1996 



for details. Figure 



[7] illustrates the Voronoi, Delaunay and contact networks for a portion of a packing of 
congruent circular disks. 

Some of the infinite packings that we will be considering in this review can only 
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FIG. 7 Geometric characterization of packings via bond networks. Left panel: Illustrations of the Voronoi 
and Delaunay tessellations for a portion of a packing of congruent circular disks. The blue and red lines 
are edges of the Voronoi and Delaunay cells, respectively. Right panel: The corresponding contact network 
shown as green lines. 

be characterized spatially via statistical correlation functions. For simplicity, consider a 
nonoverlapping configuration of N identical ^-dimensional spheres centered at the positions 
= {ri, T2, • • • , rjy} in a region of volume V in d- dimensional Euclidean space lR d . Ulti- 
mately, we will pass to the thermodynamic limit, i.e., N —> oo, V — > oo such that the number 
density p = N/V is a fixed positive constant. For statistically homogeneous sphere packings 
in W d , the quantity p n g n (r n ) is proportional to the probability de nsity for simultaneously 



findin g n sphere centers at locations r™ = {ri, r 2 , . . . , r n } in M. d ([Hansen and McDonald! . 



19861 ). With this convention, each n-particle correlation function g n approaches unity when 
all particle positions become widely separated from one another. Statistical homogeneity 
implies that g n is translationally invariant and therefore only depends on the relative dis- 
placements of the positions with respect to some arbitrarily chosen origin of the system, 
i.e., 

9n = <?n( r 12, r 13) • • • , r ln), (5) 

where = — r,j. 

The pair correlation function #2(1") is the one of primary interest in this review. If the 
system is also rotationally invariant (statistically isotropic), then <?2 depends on the radial 
distance r = |r| only, i.e., (72 (r) = gvix)- It is important to introduce the total correlation 
function h(r) = q ?(r) — 1, which, for a disorder ed packing, decays to zero for large |r| 



sufficiently rapidly (ITorquato and Stillinger 



2006bl ). We define the structure factor 5(k) for 
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a statistically homogeneous packing via the relation 



S(k) = l + ph(k), 



(6) 



where h(k) is the Fourier transform of the total correlation function h(r) = ^(r) — 1 
and k is the wave vector. Since the structure factor is the Fourier transform of an au- 



tocovariance function (involvin g; the "microscopic" density) ([Hansen and McDonald 



1986 



Torquato and Stillingerl . l2Q06bh . then it follows it is a nonnegative quantity for all k, i.e., 



S{k) > for all k. 



(7) 



The nonnegativity condition follows physically from the fact that S(k) is proportional 
to the intensity of the sca ttering of incident radiation on a many-particle system 
( IHansen and McDonald! . Il986l ). The structure factor S(k) provides a measure of the density 
fluctuations in the packing at a particular wave vector k. 



III. LESSONS FROM DISORDERED JAMMED PACKINGS OF SPHERES 



The classical statistical mechanics of hard-sphere systems has generate d a h uge col 



Boltzmann 



1898 



That 



lection of scientific publications, stretching back at least to 
collection includes examinations of equilibrium, transport, and jammed packing phe- 
nomena. With respect to the last of t hese, t h e con c ept o f a unique "random 



close packing" (RCP) state, pioneered by iBernal . 



1960 



1965 



to model the struc- 



ture of liqui d s, ha s been one of the more persistent th e mes with a venerable h istor 



(lAnonvmous 



Jullien et al 



1972 



1997 



Berrvman 



1983 



Kamien and Liu 



Gotoh and Finney 



2007 



Pouliauen et al. 



1974 



1997 



Jo drev and Tory 



Scott and Kilgour 




Tobochnik and Chapin 



1988 



Visscher and Bolsterli 



1972 



Zinchenko 



199J). Until about a 



decade ago, the prevailing notion of the RCP state was that it is the maximum density that a 
large, random collection of congruent (identical) spheres can attain and that this density is a 
well-defined quantity. This traditional view has been summarized as follows: "Ball bearings 
and similar objects have been shaken, settled in oil, stuck with paint, kneaded inside rubber 
ballo o ns-and all with no bet ter result than (a packing fraction of) . . . 0.636" ([Anonymous! . 
19721 ). iTorquato et all |2000| have argued that this RCP-state concept is actually ill-defined 



and thus should be abandoned in favor of a more precise alternative. 
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It is instructive to review briefly these developments because they will point to the need 



for a geometric-structure approach gen erally to unders t and 



ordered or not. It has been observed (ITorquato et all 120001 ) that there has existed ample 



ammed packings, whether dis- 



evidence in the literature, in the form of actual and computer-simulation experiments, to 
suggest strongly that the RCP state is indeed ill-defined and, in particular, dependent on 
the protocol used to produce the pac kings and on other system characteristics. In a classic 



experiment, 



Scott and Kilgour . 



1969 



obtained the "RCP" packing fraction value <p ~ 0.637 
by pouring ball bearings into a large container, vertically vibrating the system for sufficiently 
long times to achieve a putative maximum densification, and extrapolating the measured 
volume fractions to eliminate finite-size effects. Important dynamical parameters for this 
kind of experiment include the pouring rate as well as the amplitude, frequency and direc- 
tion of the vibrations. The shape, smoothness and rigidity of the container boundary are 
other crucial characteristics. For example, containers with curved or flat boundaries could 
frustrate or induce crystallization, respectively, in the packings, and hence the choice of 
container shape can limit the portion of configuration space that can be sampled. The key 
interactions are interparticle forces, including (ideally) repulsive hard-sphere interactions, 
friction between the particles (which inhibits densification), and gravity. The final packing 
fraction will inevitably be sensitive to these system characteristics. Indeed, one can achieve 
denser (partially and imperfectly crystalline) packings when th e particles are poured at low 
rates into horizontally shaken containers with flat boundaries ( iPouliquen et all 119971 ) . 

It is tempting to compare experimentally observed statistics of so-called RCP configu- 
rations (packing fraction, correlation functions, Voronoi statistics) to those generated on a 
computer. One must be careful in making such comparisons since it is difficult to simulate 
the features of real systems, such as the method of preparation and system characteristics 
(shaking, friction, gravity, etc.). Nonetheless, computer algorithms are valuable because 
they can be used to generate and study idealized random packings, but the final states 
are cle arly protocol- dependent. For example , a po pular rate-dependent densification algo- 



rithm ( iJodrey and To 



1985 



Jullien et al 



19971) achieves <fi between 0.642 and 0.649, a 



Monte Carlo scheme f lTobochnik and Chapinl. Il988f) g i ves <£> ~ 0.68, a differential-equation 



densific ation scheme produces 



cedure ( jVisscher and Bolsterhl . Il972[ ) yields 



0.64 ( 



Zinchenko 



19941 ). and a "drop and roll" pro- 



0.60; and each of these protocols yields 



different sphere contact statistics. 
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As noted earlier, it has been argued that these variabilities of RCP arise because it is 
an ambiguous concept, explaining why there i s no rigorous prediction o f the RCP density, 



in sp ite of attempts to estimate it (IBerrymanl . 



1983 



Gotoh and Finneyl . 



1974; 



Song et al. 



20081 ). The phrase "close packed" implies that the spheres are in contact with one another 



with the highest possible average contact number Z. This would be consistent with the 
aforementioned traditional view that RCP presents the highest density that a random pack- 
ing of close-packed spheres can possess. However, the terms "random" and "close packed" 
are at odds with one another. Increasing the degree of coordination (nearest-neighbor con- 
tacts), and thus the bulk system density, comes at the expense of disorder. The precise 
proportion of each of these competing effects is arbitrary, and therein lies a fundamental 
problem. Moreover, since "randomness" of selected jammed packings has never been quan- 
tified, the pr oportion of these com peting; effects could not be specified. To remedy these 



serious flaws, 



Torquato et all , |2000| replaced the notion of "close packing" with "jamming" 



categories (defined precisely in Sec. HV|) . which requires that each particle of a particular 
packing has a minimal number of properly arranged contacting particles. Furthermore, they 
introduced the notion of an "order metric" to quantify the degree of order (or disorder) of 
a single packing configuration. 



Using the Lubachevsky-Stillinger (LS^ 



generate jammed packings, it was shown (ITorquato et al 



(1990) molecu l ar dy namics growth algorithm to 



20001 ) that fastest particle growth 
rates generated the most disordered sphere (MRJ) packings (with <ft « 0.64; see the left 
panel of Fig. EJ), but that by slowing the growth rates larger packing fractions could be 
continuously achieved up to the densest value 7r/\/l8 ~ 0.74048. . . such that the degree 
of order increased monotonically with 0. Those results demonstrated that the notion of 
RCP as the highest possible density that a random sphere packing can attain is ill-defined, 
since one can achieve packings with arbitrarily sma ll increases in density at the expense 



of correspondingly small increases in order. This led iTorquato et al. 



2000 



to supplant the 



concept of RCP with the maximally random jammed (MRJ) state, which is defined to be 
that jammed state with a minimal value of an order metric (see Sec. IV|) . This work pointed 
the way toward a quantitative means of characterizing all packings, namely, the geometric- 
structure approach. 

We note that the same LS packing protocol that leads to a uniformly disordered jammed 
state in three dimensions typically yields a highly crystalline "collectively" jammed packing 
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FIG. 8 Typical protocols used to generate disordered sphere packings in three dimensions produce highly 
crystalline packings in two dimensions. Left panel: A three-dimensional MRJ-like configuration of 500 
spheres with <j> s a 0.64 produced using the Lubachevsky-Stillinger (LS) algorithm with a fast expansion rate 



(jTorquato et al 



2000|). Right panel: A crystalline collectively jammed configur ation (Sec. [ 



disks with <fi sa 0.88 produced using the LS algorithm with a fast expansion rate (jDonev et al 



VAl) of 1000 



2004dY 



in two dimensions. Figure El illustrates the vivid visual difference between the textures 
produced in three and in two dimensions (see Section VII for further remarks). The low- 
concentration occurrence of crystal defects in the latter is evidence for the notion that there 
are far fewer collectively jammed states for N hard disks in two dimensions compared to 
TV hard spheres in three dimensions. This distinction can be placed in a wider context by 
recalling that there is only one type of jammed state for hard rods in one dimension, and it 
is a defect-free perfect one-dimensional crystal. These cases for d = 1,2, and 3, numerical 
results for MRJ packing for d = 4, 5, and 6, and theoretical results ( Torquato and Stillingerl . 
2006bl ). indicating that packings in large dimensions are highly degenerate, suggest that the 
number of distinct collectively jammed packings (defined in Sec. HV.Aj) for a fixed large 
number N of identical hard spheres rises monotonically with Euclidean dimension d. The 
questions and issues raised by these differences in the degree of disorder across dimensions 
emphasizes the need for a geometric-structure approach, to be elaborated in the following. 
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IV. JAMMING CATEGORIES, ISOSTATICITY, AND POLYTOPES 



A. Jamming Categories 



In much of the ensuing discussion, we will treat packings of frictionless, congruent spheres 
of diameter D in R d in the absence of gravity, i.e., the "Ising model" of jammed sphere 
packings. Packing spheres is inherently a geometrical problem due to exclusion-volume 
effects. Indeed, the singular nature of the hard-sphere pair potential (plus infinity or zero 
for r < D or r > D, respectively, where r is the pair separation) is crucial because it enables 
one to be precise about the concept of jamming. Analyzing this model directly is clearly 
preferable to methods that begin with particle systems having "soft" intera ctions, which 



are th en intended to mimic packings upon passing to the hard-sphere limit (IDonev et al. 
2007dV 



Three broad and mathematically precise "jamming" categories of sphere packings can be 



distinguish ed depending on the nature of their mechanical stability (ITorquato and Stillinger 



2001 



20031 ). In order of increasing stringency (stability), for a finite sphere packing, these 



are the following: (1) Local jamming: Each particle in the packing is locally trapped by its 
neighbors (at least d+1 contacting particles, not all in the same hemisphere), i.e., it cannot be 
translated while fixing the positions of all other particles; (2) Collective jamming: Any locally 
jammed configuration is collectively jammed if no subset of particles can simultaneously be 
displaced so that its members move out of contact with one another and with the remainder 
set; and (3) Strict jamming: Any collectively jammed configuration that disallows all uniform 
volume- nonincreasing strains of the system boundary is strictly jammed. 

We stress tha t these hierarch i cal jamming; categories do n ot exhaust the universe of poss i- 



ble distinctions n Connelly et al 



1998 



Donev et al. 



2004c 



d 



Torquato and Stillinger 



20011 ). 



but they span a reasonable spectrum of possibilities. Importantly, the jamming category 
of a given sphere configuration dep ends on the boundary conditio ns employed. For exam- 
ple, hard-wall boundary conditions ( ITorquato and Stillingerl. 120011) generally y ield different 
jamming classifications from periodic boundary conditions ( IDonev et aUl2004d ). These jam- 
ming categories, which are closely r elated to the concepts of "rigid" and "stable" packings 
found in the mathematics literature ( jConnelly et all Il998l ) , mean that there can be no "rat- 
tlers" (i.e., movable but caged particles) in the packing. Nevertheless, it is the significant 
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majority of spheres that compose the underlying jammed network that confers rigidity to the 
packing, and in any case, the rattlers could be removed (in computer simulations) without 
disrupting the jammed remainder. Figure [HI shows examples of ordered locally and collec- 
tively jammed packings of disks in two dimensions within hard-wall containers. Observe the 
square-lattice packing with square hard-wall boundary conditions can only be collectively 
jammed (not strictly jammed) even in the infinite- volume limit. This is to be contrasted 
with MR J packings, where the distinction between collective and strict jamming vanishes 
in the infinite-volume limit, as discussed in Sec. IIV.BI Note that the configurations shown 
in Fig. El generated under periodic boundary conditions, are at least collectively jammed. 




FIG. 9 Illustrations of jamming categories. Leftmost panel: Honeycomb-lattice packing within a rectan- 
gular hard- wall container is locally jammed, but is not collectively jammed (e.g., a collective rotation of a 
hexagonal particle cluster, as shown, will unjam the packing). Middle panel: Square-lattice packing within 
a square hard-wall container is collectively jammed. Rightmost panel: The square-lattice packing shown 
in the middle panel can be sheared and hence is not strictly jammed. Thus, we see that the square-lattice 
packing with square hard- wall boundary conditions can only be collectively jammed even in the infinite- 
volume limit. Thus, the distinction between collective and strict jamming for such packings remains in the 
thermodynamic limit. 

To emphasize the fact that the jamming category depends on the boundary conditions 
of a packing, we tabulate whether common periodic structures are locally, collectively or 
strictly jammed. Table [I] gives the jamming classification for such packings with hard-wall 
boundary conditions. These results are compared to corresponding jamming categories for 
periodic boundary conditions in Table [Til The latter results will depend on the choice of the 
number of particles N within the fundamental cell. 
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TABLE I Classification of some of the common jammed periodic (crystal) packings of identical spheres 
in two and three dimensions, where Z denotes the contact n umber per particle and <ft is the packing 

). Here hard boundaries are applicable: in 



fraction for the infinite packing (jTorquato and Stillingerl . 



2001 



two dimensions we use commensurate rectangular boundaries and in three dimensions we use a cubical 
boundary, with the exception of the hexagonal close-packed crystal in which the natural choice is a 
hexagonal prism. 



Periodic (crystal) structures 


Locally jammed 


Collectively jammed 


Strictly jammed 


Honeycomb (Z = 3, (ft = 0.605 . . .) 


yes 


no 


no 


Kagome (Z = 4, (ft = 0.680 . . .) 


no a 


no a 


no a 


Square (Z = 4, = 0.785...) 


yes 


yes 


no 


Triangular (Z = 6, (ft = 0.907 . . .) 


yes 


yes 


yes 


Diamond (Z = 4, (ft = 0.340 . . .) 


yes 


no 


no 


Simple cubic (Z = 6, (ft = 0.524 . . .) 


yes 


yes 


no 


Body-centered cubic (Z = 8, (ft = 0.680 . . .) 


yes 


yes 


no 


Face-centered cubic (Z = 12, (ft = 0.740 . . .) 


yes 


yes 


yes 


Hexagonal close-packed (Z = 12, (ft = 0.740 . . .) 


yes 


yes 


yes 



With appropriately placed regular triangular- or hexagonal-shaped boundaries, the Kagome structure 
is locally, collectively and strictly jammed. 
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TABLE II Classification of the same periodic spher e packings described in Table UJ but for 



periodic boundary conditions (jDonev et al 



2004d ) . We chose the smallest fundamental 



cells with ./V particles for which an unjamming motion exists , if the re is one. The lattice 



vectors for each packing are the standard ones (jDonev et al. 



2004d ) 



Periodic (crystal) structures 


N 


Locally jammed 


Collectively jammed 


Strictly jammed 


Honeycomb 


4 


yes 


no 


no 


Kagome 


3 


yes 


no 


no 


Square 


2 


yes 


no 


no 


Triangular 


1 


yes 


yes 


yes 


Diamond 


4 


yes 


no 


no 


Simple cubic 


2 


yes 


no 


no 


Body-centered cubic 


2 


yes 


no 


no 


Face-centered cubic 


1 


yes 


yes 


yes 


Hexagonal close-packed 


2 


yes 


yes 


yes 



Rigorous and efficient linear-programming algorithms have been d evised to assess whet her 

2004 J dL 



a particular sphere packing is locally, collectively, or strictly jammed (jDonev et al. 



It is noteworthy that the jamming categories can now be ascertained in real-system exper- 



iments using imaging technique s that ena 



of a packing [e.g., tomography flAste et al 



and magnetic resonance imaging (IMan et al. 



e one to determine configur ational coordinates 



20061b c onfocal microscopy ( IBrujic et all 120031 ) 



20051 )]. Given these coordinates (with high 



precision), one can rigorously test the "jamming" category of the experimentally generated 
packing using the aforementioned linear programming techniques. 



B. Isostaticity 

A packing of N hard spheres of diameter D in a jammed framework in d- dimensional 
Euclidean space is characterized by the (Nd)- dimensional configuration vector of centroid 
positions R = r N = {r x , . . . , r N }. Assume that a configuration Rj represents a collectively 
jammed framework of a packing (i.e., excluding rattlers) with packing fraction 0j, where 
there are M interparticle contacts. 
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Isostatic packings are jammed packings that possess the minimal number of contacts for a 
jamming category; namely, under periodic boundary conditions, for collective jamming, M = 
27V — 1 and 3N — 2 for d = 2 and d = 3, r espectively, and for s trict jamming, M = 2N + 1 and 



3iV+3 for d = 2 and d — 3, respectively (IDonev et al. 



2005dl ). Thus, we see that the relative 



differences between isostatic collective and strict jammed packings diminish as N becomes 
large, and since the number of degrees of freedo m is essential l y equa l to Nd (depending on 
the jamming category and boundary conditions f Donev et all . 2005d )). an isostatic packing 
has a mean contact number per particle, Z, equal to 2d. Collectively /strictly jammed MRJ 
packings in the infinite-volume limit are isostatic. Note that packings in which Z = 2d are 
not necessarily collectively or strictly jammed; for example, we see from Table HT1 that the 
square and Kagome lattices (with Z = 4) and the simple-cubic lattice (with Z — 6) are 



field, and has been closely linked to "generic" or "random" packings (IMoukarzell. 



fact, as we will show, isostatic packings can be perfectly ordered; see also 



( Moukarzel . 


1995 


). In 


also 


Mari et al. 




2008. 



Packings having more contacts than isostatic ones are hyperstatic and those having fewer 



contacts than isostatic packings are hypostatic; for sphere pac kings, these 



cannot be collectively or strictly jammed in the above sense fjDonev et al 



atter pa ckings 



2007ah . The 



ter ms over constrained and underconstrained (or hypoconstrained), respectively, were used 



by 



Donev et al 



2007al to describe such packings. 



C. Polytope Picture of Configuration Space 

Full understanding of the many-body properties of ./V particles in a <i-dimensional con- 
tainer of content (volume) V is facilitated by viewing the system in its cL/V-dimensional 
configuration space. This is an especially useful approach for hard-particle models, and 
helps to understand the full range of issues concerning the approach to a jammed state. 
For the moment, we restrict attention specifically the cases of <i-dimensional hard spheres. 
When container content V is very large for fixed N, i.e., when packing fraction 0^0, the 
hard spheres are free to move virtually independently. Consequently, the measure (content) 
C of the available multidimensional configuration space is simply C ~ V . But decreasing 
V enhances the chance for sphere collisions and correspondingly reduces C, the remaining 
fraction of configuration space that is free of sphere overlaps. The amount of reduction is 
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related exponentially to the excess entropy (N, V) for the iV-sphere system: 

C(N,V) « V N exp[S^(N, V)/k B ], 



(8) 



where k B is Boltzmann's constant. 

In the low-density regime, the excess entropy admits of a power series expansion in 
covering fraction 



S&(N, V) 



Nk B 4- f \n + 1 / v ui 

n>l 

Here t>i is the volume of a particle, as indicated earlier in Eq. (J2J). The /3„ are the zrre- 
ducible Mayer cluster integral sums for n + 1 particles that determine the virial coefficient 
of order n + 1 flMayer and Mayerl . Il940l ). For hard spheres in dimensions 1 < d < 8, these 
coefficients for low orders 1 < n < 3 are known exac tly, and accurate numerical estimates 



20061 ). This power series represents a func- 



are available for n + 1 < 10 flClisby and McCoyl . 
tion of 4> obtainable by analytic continuation along the positive real axis to represent the 
thermodynamic behavior for the fluid phase from = up t o the freezing tran sition, which 
occurs at <fi ~ 0.4911 for hard spheres in three dimensions flNoya et all 120081 ). This value 



is slightly below the minimum density 



0.4937 at which col 



hard spheres is suspected first to occur (jTorquato and Stillinger 



ective jamming of d = 3 



20071 ). Consequently, the 



available configuration space measured by C(N, V) remains connected in this density range, 
i.e., any nonoverlap configuration of the iV spheres can be connected to any other one by a 
continuous displacement of the spheres that does not violate the nonoverlap condition. 
A general argument has been advanced that thermodynamic functions must ex perience a 



subt 



1964 



e , but distinctive, essential s ingularity at first-order phase transition points (jAndreevl . 



Fisher and Felderhol 



19701 ). In particular, this applies to the hard-sphere freezing 



transition, and implies that attempts to analytically continue fluid behavior into a metastable 
over-compressed state are dubious. Aside from any other arguments that might be brought 
to bear, this indicates that such extrapolations are fundamentally incapable of identifying 
unique random jammed states of the hard-sphere system. Nevertheless, it is clear that 
increasing beyond its value at the thermodynamic freezing point soon initiates partial 
fragmentation of the previously connected nonoverlap configuration space in finite systems. 
That is, locally disconnected portions are shed, each to become an individual jammed state 
displaying its own geometric characteristics. (The jammed tunneled crystals mentioned in 
Sec. IVII are examples of such localized regions near the freezing point.) We shall elaborate 
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on this point within this subsection after discussing the polytope picture of configuration 
space near jamming points. 

Consider decreasing the packing fraction slightly in a sphere packing that is at least 
collectively jammed by reducing the particle diameter by AD, 6 = AD/D <C 1, so that 
the packing fraction is lowered to <p = <pj (1 — 5) d . We call S the jamming gap or distance 
to jamming. It can be shown that there is a sufficiently small S that does not destroy the 



jamming confinement property, in the sense that the confi guration 



point R = Rj + AR 



remains trapped in a small neighborhood ^Tar around Rj (j Connelly 



19821 ). Indeed, there 



exists a range of positive values of S that depends on N and the particle arrangements that 
maintains the jamming confinement property. Let us call the threshold value at which 
jamming is lost. How does <5* scale with N for a particular dl An elementary analysis based 
on the idea that in order for a neighbor pair (or some larger local group) of particles to 
change places, the surrounding N — 2 (or N — 3, . . .) particles must be radially displaced 
and compressed outward so as to concentrate the requisite free volume around that local 
interchangeable group concludes that <5* ~ CN~ l l d , where the constant C depends on the 
dimension d and the original jammed particle configuration. 

It is noteworthy that for fixed N and sufficiently small 5, it can be shown that asymptot- 
ically (through first order in 8) the set of displacements that are accessible to the packing 



approaches a convex limiting polytope (a closed po 



( ISalsburg and Wood 



1962 



Stillinger and Salsbure 



yhedr on in high dimension) Par Q i7ar 



1969). This po lytope Par is determined 



from the linearized impenetrability equations (IDonev et al. 



2004c 



dj) and, for a fixed system 



center of mass, is necessarily bounded for a jammed configuration. This implies that the 
number of interparticle contacts M is at least one larger than the dimensionality dcs of 
the relevant configuration space. Examples of such low- dimensional polytopes for a single 
locally jammed disk are shown in Fig. [TUJ 



Importantly, for an isostatic contact network, Par is a simplex (IDonev et all l2005dl ). A 
d- dimensional simplex in M d is a closed convex polytope whose d + 1 vertices (0-dimensional 
points) do not all lie in a (d — l)-dimensional flat sub-space or, alternatively, it is a finite 
region of M d enclosed by d+1 hyperplanes [(d— l)-dimensional "faces"] (e.g., a triangle for d = 
2, a tetrahedron for d = 3 or a pentatope for d — 4). For overconstrained jammed packings 
(e.g., ordered maximally dense states), the limiting high-dimens ional polytopes have more 



faces than simplices do and can be geometrically very complex (ISalsburg and Wood! 



1962 
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FIG. 10 The polytope of allowed displacements 'Par for a locally jammed disk (light shade) trapped among 
three (left) or six (right, as in the triangular lattice) fixed disks (blue). The exclusion disks (dashed lines) 
of diameter twice the disk diameter are drawn around each of the fixed disks, along with their tangents (red 
lines) and the polytope Par they bound (black). The polytope for the isostatic (left) and overconstrained 
(right) case is a triangle and hexagon, respectively. 



Stillinger and Salsburgi Il969l ). The fact that Par is a simplex for an isostatic packing 
enables one to derive rigorous results, as we will now describe. 

Consider adding thermal kinetic energy to a nearly jammed sphere packing in the absence 
of rattlers. While the system will not be globally ergodic over the full system configuration 
space and thus not in thermodynamic equilibrium, one can still define a macroscopic pres- 
sure p for the trapped but locally ergodic system by considering time averages as the system 
executes a tightly confined motion around the particular configuration Rj. The probability 
distribution Pf(f) of the time- averaged interparticle forces / has been rigorously linked to 
the contact value r = D of the pair correlation function ^(^) defined in Section II, and this 
i n turn can be related to the distribution of simplex face hyperareas for the limiting polytope 



( iDonev et aUl2005dl ). Moreover, since the available (free) configuration volume scales in a 



predictable way with the jamming gap 5, one can s how that the reduced pressure is asymptot 



ically given by the free-volume equ ation of state ( IDonev et al 



2005d 



1962 



Stillinger and Salsburgi . Il969h 



Salsburg and Wood . 



V 



1 

5 



d 



P k B T 5 i-(H4>j? (10) 

where T is the absolute temperature and p is the number density. So far as the limiting 
polytope picture is concerned, the extremely narrow connecting filaments that in principle 
connect the jamming neighborhoods have so little measure that they do not overturn the 
free-volume leading behavior of the pressure, even as system size is allowed to go to infinity. 
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Although there is no rigorous proof yet for this claim, all numerical evidence strongly suggests 
that it is correct. Relation (fit)]) is remarkable, since it enables one to determine accurately 
the true jamming density of a given packing, even if the actual jamming point has not quite 
yet been reached, just by measuring the pressure and extrapolating to p — +00. 




Packing Fraction 



FIG. 11 The isothermal pha se behavior of th ree-dimensional hard-sphere model in the pressure-packing 
fraction plane, adapted from 



Torquato 



20021 Increasing the density plays the same role as decreasing 
temperature of a molecular liquid; see Fig. [3] Three different isothermal densification paths by which 
a hard-sphere liquid may jam are shown. An infinitesimal compression rate of the liquid traces out the 
thermodynamic equilibrium path (shown in green), including a discontinuity resulting from the first-order 
freezing transition to a crystal branch. Rapid compressions of the liquid while suppressing some degree 
of local order (blue curves) can avoid crystal nucleation (on short time scales) and produce a range of 
amorphous metastable extensions of the liquid branch that jam only at the their density maxima. 



This free- volume form has been used to estimate the equation of state along "metastable" 
extensions of the ha rd-sphere fluid up to th e infinite-pressure endpoint, assumed to be ran- 



dom jammed states (ITorquato 



1995b 



20021 ). To understand this further, it is useful to recall 



the hard-sphere phase behavior in three dimensions; see Fig. [TTJ For densities between zero 
and the "freezing" point (0 ~ 0.49), the thermodynamically stable phase is a liquid. Increas- 
ing the density beyond the freezing point results in a first-order phase transition to a crystal 
branch that begins at the melting point (</> ~ 0.55) and whose ending point is the maximally 
dense fee packing (0 m 0.74), which is a jammed packing in which each particle contacts 12 
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others f Mau and Huse 



19991 ) . However, compressing a hard-sphere liquid rapidly, under the 



constraint that significant crystal nucleation is suppressed, can produ ce a range of metastable 
branches whose density end points are random "jammed" packings (IRintoul and Torquatol . 



1996b; 



Torquato . 



20021 ). which can be regarded to be glasses. A rapid compression leads to 



a lower random jammed density than that for a slow compre ssion. The mos t r apid com- 
pressio n presumably leads to the MRJ state with ps 0.64 (jTorquatol . 120021 ). iTorquatol . 



1995aU bl reasoned that the functional form of the pressure of the stable liquid branch (which 
appears to be dominated by an unphysical pole at <fi = 1) must be fundamentally different 
from the free- volume form (ITUj) that applies near jammed states, implying that the equation 
of state is nonanalytic at the freezing point and proposed the following expression along any 
constrained metastable branch: 



V 



1 + 40g F 



1 - 4> f /4>j 



for 



< 6 < 



(11) 



pk B T TU " 1 - V i VJ 

where <pF ~ 0.491 is the packing fraction at the freezing point, gp « 5.72 is the corresponding 
value of the pair correlation function at contact, and 4>j is the jamm ing density, whose value 
will depend on which metastable path is chosen. Torquatol . ll995aU bl actually considered the 
more general problem of nearest-neighbor statistics of hard-sphere systems, which required 
an expression for the equation of state.] Unfortunately, there is no unique metastable branch 
(see Fig. [UJ because it depends on the particular constraints used to generate the metastable 
states or, in other words, the protocol employed, which again emphasizes one of the themes 
of this review. Moreover, in practice, metastable states of i dentica l spheres in M 3 have 



an inevitable tendency to crystallize (IRintoul and Torquato . 



1996bf ). but even in binary 



mixtures of hard spheres chosen to avoid crystallization the d ispersion of results and ultimate 



Kamien and Liu . 



2007 assumed 



nonuniqueness of the jammed states still apply. We note that 
the same free-volume form to fit the pressure of "metastable" states for monodisperse hard 
spheres as obtained from both numerical and experimental data to determine <pj. Their best 
fit yielded <pj = 0.6465. 

We note that density of s tates (vibrational modes) in packings o f soft spheres has been 



the subject of recent interest ( iSilbert et al. 



2005 



Wyart et al. 



20051 ) . Collective jamming in 



hard-sphere packings corresponds to having no "soft modes" in soft-sphere systems, i.e., no 
unconstrained local or global particle translations are allowed, except those corresponding 
to rattlers. Observe that it immediately follows that if a hard-sphere packing is collectively 
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jammed to first order in S, a corresponding configuration of purely soft repelling particles 
will possess quadratic modes in the vibrational energy spectrum for such a system of soft 
spheres. 



V. ORDER METRICS 



The enumeration and classification of both ordered and disordered jammed sphere pack- 
ings for the various jamming categories is an outstanding problem. Since the difficulty of the 
complete enumeration of jammed packing configurations rises exponentially with the number 
of particles, it is desirable to devise a small set of intensive parameters that can characterize 
packings well. One obvious property of a sphere packing is the packing fraction 0. Another 
important characteristic of a packing is some measure of its "randomness" or degree of dis- 
order. We have stressed that one ambiguity of the old RCP concept was that "randomness" 
was never quantified. To do so is a nontrivial challenge, but even the tentative solutions 
that have been p ut forth during the la s t decade have been profi t able not only to chara c terize 



sphere packings (jKansal et al. 



Tins. 



sett et al 



2001 



2000) 



2002b 



Torquato and Stillinger 



j ut also glas s es, simple liquids, and water (jErrington and Debenedettil . 



Errington et al. 



2003 



Torquato et al 



200C 



2002 



2003 



Truskett et al. 



2000|). 



One might argue that the maximum of an appropriate "entropic" metric would be a 
potentially useful way to characte rize the randomness of a packing and therefore the MR J 
state. However, as pointed out by 



Kansal et al. 



2002b 



a substantial hurdle to overcome in 



implementing such an order metric is the necessity to generate all possible jammed states 
or, at least, a representative sample of such states in an unbiased fashion using a "universal" 
protocol in the large-system limit, which is an intractable problem. Even if such a universal 
protocol could be developed, however, the issue of what weights to assign the resulting config- 
urations remains. Moreover, there are other fundamental problems with entropic measures, 
as we will discuss in Sec. IVIIIt including its significance for two-dimensional monodisperse 
hard-disk packings as well as polydisperse hard-disk packings with a sufficiently narrow size 
distribution. It is for this reason that we seek to devise order metrics that can be applied 
to single jammed configurations, as prescribed by the geometric-structure point of view. 

A many-body system of N particles is completely characterized statistically by its iV-body 
probability density function P(R; t) that is associated with finding the iV-particle system 
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with configuration R at some time t. Such complete information is virtually never available 
for large N and, in practice, one must settle for reduced information, such as a scalar order 
metric ip. Any order metric ip conventionally possesses the following three properties: (1) 
it is a well-defined scalar function of a configuration R; (2) it is subject typically to the 
normalization < ij) < 1; and, (3) for any two configurations R^ and R#, ^(R^) > ^>(R#) 
implies that configuration R^ is to be considered as more ordered than configuration R#. 
The set of order parameters that one selects is unavoidably subjective, given that there 
appears to be no single universal scalar measure of order. However, one can construct order 
metrics that lead to consistent results (e.g., common minima for jammed packings), as we 
will discuss after considering specific examples. 



A. Specific Order Metrics 



Many relevant order metrics have been de yised, but here w e briefl y describe only some 
of them. The bond-orientational order Qg ( jSteinhardt et all 119831 ) in three dimensions 
is defined in terms of the spherical harmonics Yi m (9i,(pi), where 9{ and ipi are the polar 
and azimuthal angles (relative to a fixed coordinate system) of the near-neighbor bond for 
particle pair %. A near neighbor could be defined as any sphere within a specified local 



radius (e.g., set by the first minimum of the pair c orrelation 



a sphere within a face-sharing Voronoi polyhedron ( jTorquato 



unctio n beyond contact) or by 



20021 ). The average of Q% over 



all of the near-neighbor bonds Nb provides a global measure of symmetries in many-particle 
systems. Of particular interest is the average when £ = 6, i.e., 

2\ 1/2 



Qt 




m=— 6 



1 N b 

i=i 



(12) 



since it reaches its maximum value for the perfect fee latt ice and is zero for a Poisson 
(uncorrelated) point distribution in the infinite-volume limit ( Rintoul and Torquato . fl996a ) 
(see Fig. [12] for the two-dimensional analog). It is seen that Q = QkIQ. I ^ lies in the closed 
interval [0,1] and therefore qualifies as an order metric (jTorquato et all |2000| ). 

A more local measure of bond-orientational order, Qeiocab can be obtained by evaluating 



the bond order at each sphere individually, and then averaging over all spheres (IKansal et al. 
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2002a), i.e., 



Qe, 



local 



1 N I A 6 1 n * 

1 V — V — Vy 

iv ^ I 13 ^ m ^ ( 

m=— 6 J i=l 



6m Pi; 



1/2 



(13) 



where n,- is the number of nearest-neighbors of sphere j. This i s analogous to the two 



Kansal et al. 



2000 



As 



dimensional definition of local bond-orientational order studied by 
noted in that work, such a local measure of order is more sensitive to small crystalline 
regions within a packing than is its global counterpart Qq, and thus avoids the possibility 
of "destructive" interference between differently oriented crystalline regions. 




FIG. 12 Two types of order, (a) Bond-orientational order contains information about the orientation of 
the vectors connecting neighboring particles in the packing of interest (left). If these orientations persist 
throughout the packing, as they do in a triangular lattice (right), the packing is considered to be perfectly 
bond-oricntationally ordered, (b) Translational order contains information about the relative spacing of 
particles in the packing of interest (left) relative to that of the densest packing at the same density (right). 



Using the fee lattice radial coordination structure as a reference system, one can define 
a translational order metric T based upon the mean occupation of sphere centers within 
thin concentric shells around each sphere in the packing as compared to the mean occu- 
pation of the same sh e lls in the fee lattice and the ideal gas at the same packing fraction 



( ITorquato et al. 



2000 



iTruskett et al. 



20001 ). i.e., 



T 



, ideal > 



i=l 



i=i 



fee 



idean 



(14) 
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Here rij is the average occupancy of the i th shell, with superscripts having the obvious 
meaning, and iV s heiis is the total number of shells employed. Figure [12] illustrates the two- 
dimensional analog of this order metric. 

All of the order metrics defined above were constructed to yield their maximum values of 
unity (when appropriately normalized) for the densest and most symmetrical (close-packed 
crystal) packing. Order metrics that are not based on any specific crystal structure have also 
been devised. For example, two different translational order metrics have been con structed 
that are based on functionals of the pair correlation function (jTruskett et all , 120001 ) . It has 
been suggested that l ocal density fluctuations within a "window" of a given size also can be 
a useful order metric (jTorquato and Stillingerl . 120031 1 . In particular, calculation of the local 
number variance for a variety of crystal, quasicrystal and "hyperuniform" disordered point 
patterns reveals that it provides a useful rank-order of these hyperuniform spatia l patterns 



at large length scales ( jTorquato and Stillinger 



2003 



Zachary and Torquatol . 120091 ). A hype 



runiform point pattern is one which the infinite-wavelength density fluctuations vanish or, 
equivalently, possesses a structure factor ff fk) (defined in Sec. II) that tends to zero in the 



limit k — > ( jTorquato and Stillinger 



20QJ). 



B. Characteristics of a Good Order Metric 



The specific order metrics have both strengths and weaknesses. This raises the question 
of what are the characteristics of a good order metric? There is clearly an enormous family 
of scalar functions that possess the aforementioned three generic properties of an order 
metric if}, but they may not necessarily be useful ones. It has been suggested that a good 
order metric should have the following additional properties (IKansal et all l2002bl ): (1) 
sensitivity to any type of ordering without bias toward any reference system; (2) ability to 
reflect the hierarchy of ordering between prototypical systems given by common physical 
intuition (e.g., perfect crystals with high symmetry should be highly ordered, followed by 
quasicrystals, correlated disordered packings without long-range order, and finally spatially 
uncorrelated or Poisson distributed particles); (3) capacity to detect order at any length 
scale; and (4) incorporation of both the variety of local coordination patterns as well as 
the spatial distribution of such patterns should be included. Moreover, any useful set of 
order metrics should consistently produce results that are positively correlated with one 
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another (iTorquatd . 120021 ; iTorquato et a/.U2000l ). The development of improved order metrics 
deserves continued research attention. 



VI. ORDER MAPS AND OPTIMAL PACKINGS 



The geometric-structure classification naturally emphasizes that there is a great diver- 
sity in the types of attainable jammed packings with varying magnitudes of overall order, 
density, and other intensive parameters. The notions of "order maps" in combination with 
the mathematically precise "jamming categories" enable one to view and characterize well- 
known packing states, such as the densest sphere packing (Kepler's conjecture) and max- 
imally random jammed (MRJ) packings as extremal states in the order map for a given 
jamming category. Indeed, this picture encompasses not only these special jammed states, 
but an uncountably infinite number of other packings, some of which have only recently been 
identified as physically significant, e.g., the jamming-threshold states (least dense jammed 
packings) as well as states b etween these and MR J. 

The so-called order map (ITorquato et all 120001 ) provides a useful means to classify pack- 
ings, jammed or not. It represents any attainable hard-sphere configuration as a point in 
the (p-ip plane. This two-parameter description is but a very small subset of the relevant pa- 
rameters that are necessary to fully characterize a configuration, but it nonetheless enables 
one to draw important conclu sions. For collective ja mming, a highly schematic order map 
has previously been proposed (ITorquato et all 120001 ) . 

Here we present a set of refined orde r maps for each of thr e e jamming categories i n R 3 (see 



Fig. [TBI) based both upon e arly work (IKansal et al 



most recent investigations (IDonev et al 



2004c 



2002b 



Torquato and Stillingerl . 



Torquato et all 120001) and the 



20071 ). Crucially, 



the order maps shown in Fig. [T3J are generally different across jamming categories and 
independent of the protocols used to generate hard-sphere configurations, and for present 
purposes include rattlers. In practice, one needs to use a variety of protocols to produce 



jammed configurations in order to popu 



the jammed regions shown in the Fig. (IKansal et al 



ate the interior anc 



to delineate the boundary of 



2002bl ). Moreover, the frequency of 



occurrence of a particular configuration is irrelevant insofar as the order map is concerned. 
In other words, the order map emphasizes a geometric-structure approach to packing by 
characterizing single configurations, regardless of how they were generated or their occur- 
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FIG. 13 Schematic order maps in the density-order (4>-^p) plane for the three different jamming categories 
in R 3 under periodic boundary conditions. White and blue regions contain the attainable packings, blue 
regions represent the jammed subspaces, and dark shaded regions contain no packings. The locus of points 
A- A' correspond to the lowest-density jammed packings. The locus of points B-B' correspond to the densest 
jammed packings. Points MRJ represent the maximally random jammed states, i.e., the most disordered 
states subject to the jamming constraint. It should be noted that the packings represented are not subject 
to rattler exclusion. 




FIG. 14 Three different optimal strictly jammed packings identified in the right most graph of Fig. 
[T3l Left panel: A: Z = 7. Middle panel: MRJ: Z = 6 (isostatic). Right panel: B: Z = 12. 



rence probability. In Fig. [131 the white and blue regions represent geometrically possible 
configurations, while the dark shaded regions are devoid of packings (e.g., maximally dense 
packings with very low order- metric values do not exist). Clearly, an appreciably reduced 
region of attainable packings will be occupied by jammed packings and, for any finite pack- 
ing, its size must decrease as the stringency of the jamming category increases. In the 
infinite-size limit (not depicted in Fig. [T3|) . the regions occupied by collectively and strictly 
jammed sets become identical. The following extremal points or loci in each jammed region 
are particularly interesting (Fig. [Ml) : 
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1. The locus of points A- A' corresponds to the lowest- density jammed packings. 
We denote by (f) m i n the corresponding jamming-threshold packing fraction. These 
packings are expected t o be characterized by a relatively high degree of order 

(jTorquato and Stillingerl . 120071 ). 

2. The locus of points B-B' correspond to the densest jammed packings with packing 

fraction </> max . 

3. The MRJ point represents the maximally random jammed state. Exclusion of rattlers 
from the MRJ state compromises its maximal irregularity; the corresponding displaced 
position in the order map involves a small reduction in packing fraction from <p « 0.64 

and a slight increase in order measure. 

4. More generally, any point along the boundary of the blue region is an extremal point, 

residing at the limit of attainability for the jamming category under consideration. 



A. Strict Jamming 

We will first discuss the strict-jamming order map. The densest sp here packings 



in three dimensions, which lie alon g the locus B-B' are strictly jammed (jDonev et al. 



2004c 



Torauato and Stillinger 



20011 ). implying that their shear moduli are infinitely large 



(jTorquato et all 120031 ). We take point B to correspond to the fee packing, i.e., it is the most 
ordered and symmetric densest packing. The other points along the line B-B' represent the 
stacking variants of the fee packing. All can be conveniently viewed as stacks of planar tri- 
angular arrays of spheres, within which each sphere contacts six neighbors. These triangular 
layers can be stacked on one another, fitting spheres of one layer into "pockets" formed by 
nearest-neighbor triangles in the layer below. At each such layer addition there are two 
choices of which set of pockets in the layer below are to be filled. Thus, all of the stacking 
variants can individually be encoded by an infinite binary sequence and therefore constitute 
an uncountabl y infinite number of maximally dense packings called the Barlow packings 
(jBarlowl . Il883l ). The most disordered subset of these is denoted by point B'. A rigorous 
proof that max = 7r/vT8 = 0.74048. . . has only recently appeared (jHalesl . 12005). In two 
dime nsions, the strictly jammed triangular lattice is the unique densest packing ( iFejes Tothl . 
19641 ) and so for d = 2 the line B-B' collapses to a single point B. 

The MRJ state is a well-defined minimum in an order map in that for a particular choice 
of jamming category and order metric it can be identified unambiguously. The MRJ concept 
is automatically compromised by passing either to the maximal packing density (fee and its 
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stacking variants) or the minimal possible density for strict jamming (tunneled crystals), 
thereby causing any reasonable order metric to rise on either side. This eliminates the 
possibility of a flat horizontal portion of the lower boundary of the jammed accessible region 
in the <j) — if) plane in Fig. [T3] (multiple MRJ states with different densities) and therefore 
indicates the uniqueness of the MRJ state in density for a particular order metric. Indeed, at 
least for collective and strict jamming in three dimensions, a variety of sensible order metrics 



produ ce an MRJ state with a packing fraction approximately equal to 0.64 (IKansal et al. 



2002bl ) (see Fig. [8]), close to the traditionally advocated density of the RCP state, and with 



an isostatic mean contact number Z = 6. This consistency among the different order metrics 
speaks to the utility of the order- metric concept, even if a perfect order metric has not yet 
been identified. However, the packing fraction of the MRJ state should not be confused with 



the MRJ state itsel f. It is possib 



e to hav e a rather ordered strictly jammed packing at this 



very same density ( IKansal et all l2002bl ). as indicated in Fig. [131 f° r example, a jammed 
but vacancy-diluted fee lattice packing. This is one reason why the two-parameter order 
map description of packings is not only useful, but necessary. In other words, density alone 
is far from sufficient in characterizing a jammed packing. 

The packings corresponding to the locus of points A- A' have received little attention 
until recently. Although it has not yet been rigorously established as such, a candidate 
for the lower limiting packing fraction min for strictly ja mmed packings is the subset o f 
"tunneled crystals" that contain linear arrays of vacancies ( Torquato and Stillinger . 12007 ). 
These relatively sparse structures are generated by stacking planar "honeycomb" layers one 
upon another, and they all amount to removal of one-third of the spheres from the maximally 
dense structures with packing fraction max . Consequently, min = 20 max /3 = 0.49365. . .. 
Every sphere in a tunneled crystal contacts 7 immediate neighbors in one of two possible 
coordination geometries, and all of the stacking variants exhibit some form of long-range 
order. It is appropriate to view the two families of maximum- density and of minimum- density 
strictly jammed packings as structural siblings of one another. Note that jammed packings 
can trivially be created whose densities span the entire range between these extremal cases 
simply by filling an arbitrary fraction of the vacant sites in any one of the tunneled structures. 
The dashed lines joining the points A to B and points A' to B' shown in Fig. [13] are the 
result of sequentially filling the most ordered and disordered tunneled crystals with spheres 
until the filling process ends with the most ordered and disordered densest Barlow packings, 
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respectively. Interestingly, the tunneled crystals exist at the edge of mechanical stability, 
since removal of any o ne sphere from the in t erior w ould cause the entire packing to collapse. 



It is noteworthy that 



Burnell and Sondhi 



2008 have shown that an infinite subclass of 



the tunneled crystals has an underlying topology that greatly simplifies the determination 
of their magnetic phase structure for nearest-neighbor antiferromagnetic interactions and 
O(N) spins. 

It should come as no surprise that ensemble methods that produce "most probable" 
configurations typically miss interesting extremal points in the order map, such as the locus 
of points A- A' and the rest of the jamming- region boundary, including remarkably enough 
the line B-B' . However, numerical protocols can be devised to yield unusual extremal 
jammed states, as discussed in Sec. IVII[ for example. 



Observe that irregular jammed packings c an be createc 



of packing fraction 0.64 < < 0.74048... (IKansal et al. 



in the e ntire non-trivia l range 



2002b 



Torquato et all |2000| ) 



using the LS algorithm. Thus, in the rightmost plot in Fig. [13l the MRJ-B' portion of the 
boundary of the jammed set, possessing the lowest order metric, is demonstrably achievable. 
Until recently, no algorithms have produced disordered strictl y jammed packings to the 



left of the MRJ point. A new algorithm described elsewhere (jTorquato and Jiaol . I2010d ) 
has indeed yielded such packings with <p m 0.60, which are overconstrained with Z m 6.4, 
implying that they are more ordered than the MRJ state (see Sec. I VI 1 1 for additional details). 
The existence of disordered strictly jammed packings with such anomalously low densities 
expands conventional thinking about the nature and diversity of disordered packings and 
places in a broader context those protocols that produce "typical" configurations. 

Indeed, there is no fundamental reason why the entire lower boundary of the jammed set 
between the low-density jamming threshold and MRJ point cannot also be realized. Note 
that such low-density disordered packings are not so-called "random loose" packings, which 
are even less well-defined than RCP states. For example, it is not clear that the former are 
even collectively jammed. A necessary first step would be to classify the jamming category 
of a random loose packing (RLP), which has yet to be done. Therefore, in our view, the 
current tendency in the li terature to put so- called RCP and RLP on the same footing as far 



20081 ) is premature at best. 



as jamming is concerned (jSong et all 

In IR 2 , the so-called "reinforced" Kagome packing with precisely 4 contacts per particle 
(in the infinite-packing limit) is evidently the lowest density strictly jammed subpacking of 
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the triangular lattice packing (IDonev et al. 



2004ch with 



= .68017.... Note 



that this packing has the isostatic contact number Z = 4 and yet is an ordered packing, 
which runs cou nter to the prevale nt notion that isostaticity is a consequence of "genericity" 



or randomness (IMoukarzel 



1998|). 



B. Collective and Local Jamming 



Observe that the locus of points B-B' is invariant under change of the jamming category, 
as shown in Fig. [13j This is not true of the MRJ state, which will generally have a differ- 
ent location in the local-jamming and collective-jamming order maps. Another important 
distinction is that it is possible to pack spheres subject only to the weak locally-jammec 
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criterion, so that the resulting packing fraction is arbitrarily close to zero (iBoroczky 
Stillinger et all 120031 1 . But demanding either collective jamming or strict jamming evidently 
forces <p to equal or exceed a lower limit min that is well above zero. 



C. Broader Applications to Other Condensed States of Matter 



Although methods for character i zing s t ructural order in regula r 



established (lAshcroft and Mermin 



1976 



Chaikin and Lubensky 



cryst alline solids are well 



19951 ) . similar techniques 



for noncrystalline condensed states of matter are not nearly as advanced. The notions of 
order metrics and order maps have been fruitfully extended to characterize the degree of 
structural order in condensed phases of matter in which the constituent molecules (jammed 
or not) possess both attractive and repulsive interactions. This includes the determination 



of the order map s of models of simple 



( Errii 


agton et al. 


2003 


Truskett et al. 


2001; 


Errington et al.. 


2002 


), and moc 



iquids, glasses and crysta 



2000), models of water flErrington anc 



s with isotropic interactions 



Debenedettil . 



2003|). 



VII. PROTOCOL BIAS, LOSS OF ERGODICITY, AND NONUNIQUENESS OF 
JAMMED STATES 

A dilute system of N disks or spheres is free to reconfigure largely independently of 
particle-pair nonoverlap constraints. However, as <fi increases either as a result of compression 
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or of particle size growth, those constraints consume larger and larger portions of the dN- 
dimensional configuration space, making reconfiguring more and more difficult. Indeed, 
the available subspace begins to fracture, prod ucing isolated "islands" th at each eventually 



collapse with increasing <fi into jammed states (jSalsburg and Woodl . 119621 ). This fracturing, 
or disconnection, implies dynamical non-ergodicity. Owing to permutation possibilities for 
N identical objects, each disconnected region belongs to a large family of essentially N\ 
equivalent regions. But the configuration space fracturing has even greater complexity in 
that the number of inequivalent such families and their jamming limits rises exponentially 
with N in the large-system asymptotic limit. The jamming-limit values for the families 
vary over the ranges indicated in Fig. [T3]for collective and strict jamming. 



Pre-Jammed States 




FIG. 15 A schematic of the disjoint set of nearly jammed packings that develop in multidimensional 
configuration space as covering fraction (j) increases. The two axes represent the collection of configurational 
coordinates. Note that each individual region approaches a polytope in the jamming limit, as discussed in 
Sec. UTO 

Figure [15] offers a simple schematic cartoon to illustrate this <iiV-dimensional discon- 
nection feature. Several allowed regions with different sizes and shapes are shown. Their 
boundaries consist of sets of slightly curved hypersurfaces, each of which corresponds to a 
particle pair contact, or contact with a hard wall if present. Particle growth or system com- 
pression causes hypersurfaces (numbering at least dN + 1 for hard walls) to move inward, 
reducing region content toward zero. The larger the region shown, the larger should be 
understood its jamming <fi value. The basic issue involved in either laboratory or computer 
experiments is how and why the various jamming protocols used populate the disconnected 



39 



regions. Presumably any given algorithm has associated with it a characteristic set of occu- 
pation weights, leading in turn to well-defined averages for any property of interest, including 
packing fraction and any chosen order metric ift. The fact that these averages indeed vary 
with algorithm is a major point of the present review. 

Ensemble methods have been invoked to attach special significance to so-called "typical" 
or "unique" packings because of their frequency of occurrence in the specific method em- 
ployed. In particular, significance has been attached to the so-called unique J (jammed) 
point, w hich is suggested to correspond to the onset of collective jamming in soft sphere 



systems fjO'Hern et al 



20031 ). The order maps described in Sec. |V] as well as the ensuing 



discussion demonstrate that claims of such uniqueness overlook the wide variability of pack- 
ing algorithms and the distribution of configurations that they generate. Individual packing 
protocols (numerical or experimental) produce jammed packings that are strongly concen- 
trated in isolated pockets of configuration space that are individually selected by those 
protocols. Therefore, conclusions drawn from any particular protocol are highly specific 
rather than general in our view. 

Indeed, one can create protocols that can lead to jammed packings at any preselected 
density with a high probability of occurrence anywhere over a wide density range. Unless it 
were chosen to be highly restrictive, a typical disk or sphere jamming algorithm applied to a 
large number iV of particles would be capable of producing a large number of geometrically 
distinguishable results. In particular, these distinguishable jammed configurations from 
a given algorithm would show some dispersion in their and ip values. However, upon 
comparing the distributions of obtained results for a substantial range of particle numbers 
iV (with fixed boundary conditions), one must expect a narrowing of those distributions with 
increasing iV owing to operation of a central limit theorem. Indeed, this narrowing would 
converge individually onto values that are algorithm-specific, i.e., different from one another. 
Figure [T6l provides a clear illustration of such narro wing with respect t o distributions, with 



evident variation over algorithms, as obtained by Uiao et all l2010bl The examples shown 
contrast results for two distinctly different sphere system sizes (~ 250 and ~ 2500 particles), 
and for two different algorithms that have results for disordered jammed packings converging, 
respectively, onto packing fractions of about 0.60, 0.64, 0. 68 and 0.72 The his t ogram for the 
lowest density was produced using the new algorithm (ITorquato and Jiaol . l2010d ) noted 
in Section IVI[ while the other two histograms were generated using the LS algorithm. We 
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FIG. 16 Packing protocols can be devised that lead to strictly jammed states at any specific density with 
a high probability of occurrence anywhere over a wide density range. Shown are histograms of jammed 
pa ckings that are cen tered around four different packing fractions: </> ps 0.60, 0.64, 0.68 and 0.72, as obtained 



by 



Jiao et al 



2010U The distributions become narrower as the system size increases. 



stress once again that any temptation to select a specific <f> value as uniquely significant (e.g., 
0.64) is primarily based on inadequate sampling of the full range of algorithmic richness and 
diversity that is available at least in the underlying mathematical theory of sphere jamming. 



VIII. ATTRIBUTES OF THE MAXIMALLY RANDOM JAMMED STATE 



The MRJ state u n der t he strict-jamming constraint is a prototypical glass 



(jTorquato and Stillingerl . 120071 ) in that it is maximally disordered wit hout any long 



range order and perfectly rigid (the elastic moduli are indeed unbounded (jTorquato et al. 
2003T) ) . This endow s such packings, which are isostatic {Z = 2d), with special attributes 



( IDonev et aUl2005al Jdl). For example, the pair correlation function g2{r) (which provides the 
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<TG. 17 Pair statist ics for packings in the immediate neighborhood of the three-dimensional MRJ state 
2002bJ) with <f> w 0.64. Left panel: Computatio nal data on the pair correlation function 



( Kansal et al 



2005dl ) of diameter D. It is 



.92(7) versus r/D — 1 averaged over 5 packings of 10,000 spheres (jDonev et al. 
normalized so that at large distances it tends to unity, indicating no long-range order. The split second peak, 
the discontinuity at twice the sphere diameter, and the divergence near contact are clearly visible. Right 
panel: The corresponding s tructure factor S(k) as a function of the dimcnsionless wavenumber kD /(2tt) for 



2005a)). The inset shows the linear ( |k|) nonanalytic behavior at 



a million particle packing (jDonev et al 
k = 0. 

distri bution of pair distances) of three-dimensional MRJ packings possesses a split second 
peak ( IZallenl . Il983l ). with a prominent discontinuity at twice the sphere diameter, as shown 
in the left panel of Fig. [T71 which is a well-known characteristic of disordered jammed pack- 
ings. The values r = \/3D and r = 2D are highlighted in the figure, and match the two 
observed singularities. Interesting ly, an integrable pow er-law divergence (l/(r/D— l) a with 
a ~ 0.4) exists for near contacts f Donev et all l2005dl ). No peaks are observed at r = \/2D 



or r = V5D, which are typical of crystal packings, indicating that there is no detectable 
undistorted cry stal ordering i n the packing. We note that in a computational study of stiff 
"soft" spheres (ISilbert et all 120021 ). a nearly square-root divergence for near contacts was 
found. 

The MRJ state possesses unusual spatial density fluctuations. It was conjectured 



( ITorquato and Stillingerl . 120031 ) that all strictly jammed saturated packings of congruent 



spheres (disordered or not) are hyperuniform, i.e., infinite- wavelength density fluctuations 
vanish or, equivalently, the structure factor S(k) vanishes in the limit k — > 0. (Recall that a 
saturated packing is one in which no space exists to insert additional particles.) Disordered 
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hyperuniform point distributions are uncommon. Not only was this conjecture verified nu- 
merically for an MRJ-like state using a million-particle packing of monodisperse spheres, but 
it was sho wn that the structur e factor has an unusual nonanalytic linear dependence near 
the origin (IDonev et all l2005af ). namely, S(k) ~ \k\ for k — >• 0, or equivalently, a quasi-long- 



ranged tail of the total pair correlation function h(r) 



This same linear nonanalytic 



behavior of S(k) near the origin is also found in such diverse thre e-dimensional systems 



as th e early Universe ( 



956 



Peebles 



Reatto and Chester. 



flScardicchio et al 



2009 



.9931 ). ground state of liquid helium (IFeynman and Cohenl . 



19671) and nqninter cting spin-polaried fermionic ground states 



Torquato et al. 



20081 ). The generalization of the aforementioned 



conjecture that describes when strictly jammed saturated pa ckings of noncqngruen t sphere 



Zachary et al 



2010 



Specif- 



packings as well as other particle shapes has been given by 
ically, the void spaces of general MRJ packings are highly constrained by the underlying 
contact network, which induce hyperuniformity and quasi-long-range behavior of the two- 
point probability function for the void phase. 

The quasi-long-range behavior of g2{r) as well as the aforementioned pair-correlation 
features distinguish the MRJ state strongly from that of the equilibrium hard-sphere fluid 
( [Hansen and McDonald! . Il986f ). which is characterized by a structure factor that is analytic 
at k = and thus has a pair correlatio n function tha t decays exponentially fast to unity for 



large r. Consequently, early attempts (IBernal . 
model liquid structure were imprecise. 

It shoul d be recognized that MRJ-like sphere packings created in 



19601 ) to use disordered jammed packings to 



algorithms (iDonev et al. 



2005a 



O'Hern et al. 



2002 



pract ice via computer 



Torquato et all 120001 ) or actual exper- 



iments may contain a small concentration of rattlers, the average concentration of which is 
protocol-dependent. The packings leading to the data depicted in Fig. [T7| contain between 
2-3 percent rattlers. Thus, the hyperuniformity property of the MRJ state requires that the 
rattlers be retained in the packing. 



f Jullien et al. 


1997; 


Torquato. 


2002) 



dimensional analogs of three-dimensional computational and experimental protocols that 
lead to putative RCP states result in packings of identical disks that are highly crystalline, 
forming rather large triangular coordination domains (grains). Such a 1000-particle packing 
with cf) w 0.88 is depicted in the right panel of Fig. [HJ and is only collectively jammed at 
this high density. Because such highly ordered packings are the most probable outcomes 
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for these typical protocols, "entropic measures" of disorder would identify these as the 
most disordered, a misleading conclusion. An appropriate order metric, on the other hand, 
is capable of identifying a particular configuration (not an ensemble of configurations) of 
considerably lower density (e.g., a jammed vacancy-diluted triangular lattice or its multi- 
domain variant) that is consistent with our intuitive notions of maximal disorder. However, 
typical packing protocols would almost never generate such disordered disk configurations 
because of their inherent implicit bias toward undiluted crystallization. Note that the same 
problems persist even for polydisperse disk packings provided that the size distribution is 
sufficiently narrow. 



Importantly, previous 



Donev et al 



V reported 



low packing fractions of 0.82 - 0.84 for so-called 



RCP disk arrangeme nts (iBerrymanl . Il983l ) were found not even to be collectively jammed 



2004d ). This conclusion clearly demonstrates that the distinctions between 



the different jamming categories are crucial. Moreover, the geometric-structure approach to 
jamming reveals the basic importance of collective motions potentially involving an arbitrar- 
ily large number of particles. Therefore, methods that assume collective jamming based only 
on packin g fraction and loc al criteria, such as nearest-neighbor coordination and Voronoi 
statistics (jSong et all 120081 ). are incomplete. 



IX. PACKINGS OF SPHERES WITH A SIZE DISTRIBUTION 



Polydispersity in the size of the particles constitutes a fundamental feature of the mi- 
crostructure of a wide class of dispersions of tech nological i mpor tance, including those 



i nvolved in comp osite soli d propellant combu stion (jKersteinl . 119871 ). sintering of powders 



(IRahaman 



ticula te composite materials ( 
19740 . 



Christensen 



19951 ). colloids (IRussel et q/.LIl989f). tr ansport and mechanical prop erties of par 



19791 ). and flow in packed beds (jScheideggerl . 



The spheres generally possess a distribution in radius R characterized by a probability 
density f(R) that normalizes to unity, i.e., 



f(R)dR = 1. 



The average of any function w(R) is defined by 



(w(R)) 



w{R)f{R)dR. 



(15) 



(16) 
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1.0 2.0 

R/<R> 



1.0 2.0 

R/<R> 



FIG. 18 Frequently used hard-sphere size distributions. Left panel: Schulz size distribution for 
several values of m. Right panel: Log-normal size distribution for several values of j3. 



The overall packing fraction of the system is defined as 



<j> = p{v 1 {R)), 



(17) 



where p is the total number density, Vi(R) is given by 02J and (vi(R)) is the average sphere 
volume defined by 



(vi(R)) 



7T 



d/2 



(R d ). 



(18) 



r(l + d/2) 

There is a variety of choices for the size distribution f(R) that deserve consideration. 
Two continuous probabilit y densities th at have been wi dely used to ch aracterize physical 
phenomena are the Schulz ( Schulz . 19391 ) and log- normal (Cramer, 1954 ) distributions. The 
Schulz distribution is defined as 



f(R) 



m + 1 



m+l 



R m exp 



-(m + l)R 

W) . 



(19) 



r(m + 1) V (R) 

where T(x) is the gamma function. When the parameter m is restricted to nonnegative 
integer values, T(m + 1) = ml, and the nth moment of this distribution is given by 

(m + n)l 1 



(R n ) 



(R) n - 



(20) 



mi (m + 1)" 

By increasing m, the variance decreases, i.e., the distribution becomes sharper. In the 
monodisperse limit m — > oo, f(R) 5(R — (R))- The case m = gives an exponential 
distribution in which many particles have extremely small radii. By contrast, the log-normal 
distribution is defined as 



f(R) 



R^h^'- 



; exp 



[HR/(R))f 

2/3 2 



(21) 
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where /3 2 = ((lni?) 2 ) — (lni?) 2 . The quantity InR has a normal or Gaussian distribution. 
The nth moment is given by 



(R n ) = exp(n 2 /372) (R} n . 



(22) 



As f3 2 — > 0, f(R) — > S(R — (R)). Figure [TBI shows examples of the Schulz and log-normal 
size distributions. 

One can obtain corresponding results for spheres with M discrete different sizes from the 
continuous case by letting 



M 



f{R) = Y,-A R 
i=i p 



Ri 



(23) 



where pi and Ri are number density and radius of type-z particles, respectively, and p is the 
total number density. Therefore, the overall volume fraction using ( fl7l) is given by 

M 

= ^0« (24) 



i=i 



where 



W =p i v 1 (R i ) 



(25) 



is the packing fraction of the zth component. 

Sphere packings with a size distribution exhibit intriguing structural features, some of 
which are only beginning to be understood. It is known, for example, that a relatively small 
degree of polydispersity ca n suppress the disorder- order phase transition seen in monodis- 
perse hard-sphere systems ([Henderson et all 119961 1 . Interestingly, equilibrium mixtures of 
small and large hard spheres can "phase separate" (i.e., the small and large spheres demix) at 
sufficiently high densities but the precise nature of su ch phase transi t ions h as not yet been 
established and is a subject of intense interest; see (IDijkstra et all Il999l ) and references 
therein. 

Our main interest here is in dense polydisperse packings of spheres, especially jammed 
ones. Very little is rigorously known about the characteristics of such systems. For example, 
the maximal overall packing fraction of even a binary mixture of hard spheres in M d , which 



we denote by 



;(2) 



for arbitrary values of the mole fractions and radii R\ and R^ is unknown, 



not to mention the determination of the corresponding structures. However, one can bound 
0mL from above and below in terms of the maximal packing fraction 0mL f or a mono dispers e 



sphere packing in the infinite-volume limit using the following analysis of (jTorquato 



20021 ) 
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It is clear that 0mL is bounded from below by the maximum packing fraction 0mL- The 
lower bound 0mL > 0mL is independent of the radii and corresponds to the case when 
the two components are completely phase separated (demixed), each at the packing fraction 
0max- Moreover, one can bound <j)mlx from above in terms of the monodisperse value </>mix 
for arbitrary values of R\ and i?2- Specifically, consider a wide separation of sizes (R\ <C R2) 
and imagine a sequential process in which the larger spheres are first packed at the maximum 
density 0mL for a monodisperse packing. The remaining interstitial space between the larger 
spheres can now be packed with the smaller spheres at the packing fraction 0mL provided 
that R1/R2 —> 0. The overall packing fraction in this limit is given by 1 — (1 — 0mix) 2 , which 
is an upper bound for any binary packing. Thus, 0mL < 1 — (1 — tt/v^12) 2 ~ 0.991 for d = 2 
and </>mix < 1 — (1 — 7r/vT8) 2 ~ 0.933 for d = 3, where corresponds to the maximal 
packing fraction in two or three dimensions, respectively. 

The same argume nts extend to systems of M different hard spheres with radii 
Ri, R2, ■ ■ • , Rm in M. d (jTorquatol . 120021 ). Specifically, the overall maximal packing fraction 
0max of such a general mixture in ~R d [where <p is defined by ( jTTj) with (j23l) ] is bounded from 
above and below by 

€1 < & < 1 - (1 - 0mL) M (26) 

The lower bound corresponds to the case when the M components completely demix, each at 
the density <f)mlx- The upper bound corresponds to the generalization of the aforementioned 
ideal sequential packing process for arbitrary M in which we take the limits R1/R2 — > 0, 
R2/R3 0, • • • , Rm-i/Rm 0. Specific nonsequential protocols (algorithmic or other- 
wise) that can generate structures that approach the upper bound ( 1261) for arbitrary values 
of M are currently unknown and thus the development of such protocols is an open area 
of research. We see that in the limit M — > 00, the upper bound approaches unity, cor- 
responding to space-fil 



ing polydisperse spheres with an infinitely wide separation in sizes 
( IHerrmann et aUll990l ). Furthermore, one can also imagine constructing space-filling poly- 
disper se spheres with a continuous size distribution with sizes ranging to the infinitesimally 
small (jTorquatol . 120021 ) . 

Jammed binary packings have received some attention but their characterization is far 
from complete. Here we briefly note work concerned with maximally dense binary packings 
as well as disordered jammed binary packings in two and three dimensions. Among these 
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FIG. 19 One large d isk and two sm all disks in mutual contact provide the densest local arrange- 



ment of binary disks (jFlorian 



1960). The intersection of the shaded triangle with the three disks 



ira 2 + 2(1 — a 2 )arcsin 



yields the local packing fraction <f>u 



a 



l + a 



2a(l + 2a) 



l/T 



, where a = Rs/Rl- 



cases, we know most about the determination of the maximally dense binary packings in 
IR 2 . Let Rs and Rl denote the radii of the small and large disks (R$ < Rl), the radii ratio 
a = Rs I Rl and Xs be the number fraction of small disks in the entire packing. Ideally, it is 
desired to obtain max as a function of a and xs- In pract ice, we have a sketchy understanding 



1964 



has reported a number of candidate 



of the surface denned by (j) max (a, xs)- iFejes Tothl . 
maximally dense packing arrangements for certain values of the radii ratio in the range 
a > 0.154701 .... Maximally dense binary disk packings have been also in v estiga ted to 
determine the stable crystal phase diagram of such alloys (ILikos and Henleyl . Il993l ). The 
determination of max for sufficiently small a amounts to finding the optimal arrangement of 
the small disks within a tricusp: the nonconvex cavity between three close-packed large disks. 
A particle-growth Monte Carlo algorithm was used to generate the densest arrangements 
of small identical disks (ranging in number from one through 19) within such a tricusp 



(lUche et al. 



20041 ) . All of these results can be compared to a relatively sharp upper bound 
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on 



given by 



no, 2 + 2(1 — a 2 )arcsin 



a 



*- £ *" = mTTW ' (27) 

The fraction <fiu corresponds to the densest local packin g arrangement for three binary disks 
shown in Fig. [19] and hence bounds max from above ( Florian , 196oL Inequality (127|) also 
applies to general multicomponent packings, where a is taken to be the ratio of the smallest 
disk radius to the largest disk radius. 

The most comprehensive study of the densest possible packings of binary spheres in 
R 3 as well as more general s ize-discrete mixtures has been reported in a recent paper by 



Hudson and Harrowell 



2008 



These authors generated candidate maximally dense polydis- 
perse packings based on filling the interstices in uniform three-dimensional tilings of space 
with spheres of different sizes. They were able to fin d for certain size ratios and com positions 
a number of new packings. The reader is referred to iHudson and Harrowell 120081 for details 
and some history on the three-dimensional problem. 

One of the early numerical investigations of disordered jammed pa ckings of binary disks in 
R 2 an d binary spheres in IR 3 employed a "drop and roll" procedure (IVisscher and Bolsterlil . 



19721 ). Such numerical protocols and others (jOkubo and Odagakil . 120041 ). in which there 



is a preferred direction in the system, tend to produce statistically anisotropic packings, 
which exhibit lower densities than than those gen erated by packing protocols that yield sta- 
tistically isotropic packings (IDonev et all l2004q ). It is not clear that the former packings 
are collectively jammed. In two dimensions, one must be especially careful in choosing a 
sufficiently small size ratio in order to avoid the tendency of such packings to form highly 
crystalline arrangements. The LS algorithm has been used successfully to g enerate disor 
dered strictly jammed packings of binary disks with <fi w 0.84 and or 1 = 1.4 (IDonev et al. 



20061 ). By explicitly constructing an exponential number of jammed packings of binary disks 



with densities spanning the spectrum 
separated cr ystal, it has been argued ( 



rom the accepted amorpho us glassy state to the phase- 



Donev et al. 



2006 



2007bl ) that there is no "ideal glass 



transition" fjParisi and Zamponil . 120051 ). The existence of an ideal glass transition remains 
a hotly debated topic of research. 



In three dimensions, it was shown by 



Schaertl and Sillescu 



1994! that increasing poly- 



dispersity increases the packing fraction over the monodisperse value that an amorphous 
hard-sphere system can possess. The LS algorithm has been extended to generate jammed 
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sphere packings with a polydispersity in size (IKansal et all |2002aj). It was applied to show 
that disordered packings with a wide range of packing fractions that e xceed 0.64 and varying 



degrees of disorder can be achieved; see also 



Chaudhuri et al. 



2010 



Not surprisingly, the 



determination of the maximally rando m jammed (MRJ) s tate for an arbitrary polydisperse 



sphere packing is a wide open question. 



Clusel et al. 



2009 



have carried out a beautiful series 



of experiments to understand polydisperse random packings of spheres. Specifically, they 
produced three-dimensional random packings of frictionless emulsion droplets with a high 
degree of size polydispersity, and visualize and characterize them using confocal microscopy. 

The aforementioned investigations and the many conundrums that remain serve to il- 
lustrate the richness of polydisperse jammed packings but further discussion is beyond the 
scope of this review. 

X. PACKINGS OF NONSPHERICAL PARTICLES 



Jammed packing characteri stics become co n siderably more comp l ex by allow i ng for 



nonspherica 


Darticle shapes (Betke and Hen! 


:. 2000: 


Chai 


tin et al. 


200' 


7 ; 


Chen. 


2008; 


Chen et al. 


2010: 


Conwav and Tore 


uato, 


2006; 


Donev et al. 


2004a. 2007c 


t. 


2004b. 2005blc: 


Haii-Akbari et al. 


2009: Jiao et al. 


2010a 


Kallus et al. 


20 


L0 


: Man et al 




2005; 


Roux 


2000; 


Torquato and Jiao 


. 2009a b. 2010b 


Williams and Philipse 


• 


2003). We focus here on the 



latest developments in this category, specifically for particle shapes that are continuous 
deformations of a sphere (ellipsoids and superballs) as well as polyhedra. Nonsphericity 
introduces rotational degrees of freedom not present in sphere packings, and can dramati- 
cally alter the jamming characteristics from those of sphere packings. [We note in passing 
that there has been deep and productive examinati on of the equilibrium phase behavior and 



transport properties of hard nonspherical particles ((Allen 



1993 



Frenkel et al 



1984 : 



Yatsenko and Schweizer 



20081 ).] 



Frenkel and Maguire 



198 



3; 



Very recent developments have provided organizing principles to characterize and classify 
j ammed packings o f nonsp herical particles in terms of shape symmetry of the particles 
( ITorquato and Jiaol . |2009aJ). We will elaborate on these principles in Sec JX.CI where we 
discuss packings of polyhedra. We will begin the discussion by considering developments 
within the last several years on ellipsoid packings, which has spurred much of the resurgent 
interest in dense packings of nonspherical particles. 
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A. Ellipsoid Packings 



One simple generalization of the sphere is an ellipsoid, the family of which is a contin- 
uous deformation of a sphere. A three-dimensional ellipsoid is a centrally symmetric body 
occupying the region 

© 2 +(?r+(?) 2 ^ <*> 



where (i = 1,2,3) are Cartesian coordinates and a, b and c are the semi-axes of the 
ellipsoid. Thus, we see that an ellipsoid is an affine (linear) transformation of the sphere. A 
spheroid is an ellipsoid in which two of the semi-axes are equal, say a = c, and is a prolate 
(elongated) spheroid if b > a and an oblate (flattened) spheroid if b < a. 

Figure [20] shows how prolate and oblate spheroids are obtained from a sphere by a linear 
stretch and shrinkage of the space along the axis of symmetry, respectively. This figure 
also illustrates two other basic points by inscribing the particles within the smallest circular 
cylinders. The fraction of space occupied by each of the particles within the cylinders is an 
invariant equal (due to the affine transformations) to 2/3. This might lead one to believe that 
the densest packing of ellipsoids is given by an affine transformation of one of the densest 
sphere packings, but such transformations necessarily lead to ellipsoids that all have exactly 
the same orientations. Exploiting the rotational degrees of freedom so that the ellipsoids are 
not all required to have the same orientations turns out to lead to larger packing fractions 
than that for maximally dense sphere packings. Furthermore, because the fraction of space 
remains the same in each example shown in Figure [2Ql the sometimes popular notion that 
going to the extreme "needle-like" limit (b/a — > oo) or extreme "disk-like' limit (b/a — > 0) 
can lead to packing fractions <fi approaching unity is misg uided. 



Experiments on M&M candies (spheroidal particles) (IDonev et al. 



2004a 



2005J) as well as numerical results produced by a modified LS algorithm (IDonev et al 



Man et al 



2005bU d) found MRJ-like packings with packing fractions and mean contact numbers that 
were higher than for spheres. This led to a numerical study of the packing fraction and 
mean contact number Z as a function of the semi-axes (aspect) ratios. 

The results were quite dramatic in several respects. First, it was shown that and Z, as 
a function of aspect ratio, each have a cusp (i.e., non-differentiable) minimum at the sphere 
point, and 4> versus aspect ratio posses ses a density maximu m; see Fig. [2TJ which shows the 



more refined calculations presented in fjDonev et al. 



2007al ). The existence of a cusp at the 
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(i) (ii) (iii) 

FIG. 20 A sphere inscribed within the smallest circular cylinder (i) undergoes a linear stretch and 
shrinkage of the space along the vertical direction leading to a prolate spheroid (ii) and an oblate 
spheroid (iii), respectively. This linear transformation leaves the fraction of space occupied by the 
spheroids within the cylinders unchanged from that of the fraction of the cylinder volume occupied 
by the sphere, equal to 2/3. 



0.74 



0.72 




0.68 
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IG. 21 Density (ft versus aspect ratio a for MRJ packings of 10,000 ellipsoids as obtained in 



(Donev et al 



2007al ). The semi-axes here are l,a",a. The inset shows the mean contact number 
Z as a function of a. Neither the spheroid (oblate or prolate) or general ellipsoids cases attain 
their isostatic values of Z = 10 or Z = 12, respectively. 
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sphere point runs counter to the prevailing expectation in the literature that for "generic" 
(disordered) jammed frictionless particles the total number of (independent) constraints 
equals the total number of degrees of freedom df, implying a mean contact number Z = 2df 
(df — 2 for disks, df = 3 for ellipses, df = 3 for spheres, df = 5 for spheroi ds, and df = 6 for 
general ellipsoids). This has been referred to as the isostatic conjecture ( lAlexanderl . Il998l ) 
or isocounting conjecture (IDonev et a/I l2007al ). Since df increases discontinuously with the 
introduction of rotational degrees of freedom as one makes the particles nonspherical, the 
isostatic conjecture predicts that Z should have a jump increase at aspect ratio a = 1 
to a value of Z — 12 for a general ellipsoid. Such a discontinuity was not observed by 



Donev et al. 



2004al . rather, it was observed that jammed ellipsoid packings are hypostatic, 
Z < 2df, near the sphere point, and only become nearly isostatic for large aspect ratios. 
In fact, the isostatic conjecture is only rigorously true for amorphous sphere packings after 



removal of ratt 



sub-isostatic) ( IDonev et al. 



ers; generic nonspherical-particle packings should generally be hypostatic (or 



2007a 



Roux 



2000|). 



Until recently, it was accepted that a sub-isostatic or hypostatic packi ng of nonspherica l 



particles cannot be rigid (jammed) due to the existence of "floppy" modes (lAlexanderl . 119981 ) . 
which are unjamming motions (mechanisms) derived within a linear theory of rigidity, i.e., 
a first-order analysis in the jamming gap 5 (see Sec. HV.Cj) . The observatio n that terms 



Roux 



2000 



of order higher than first generally need to be considered was emphasized by 
but this analysis was only developed for spheres. It has recently been rigorously shown that 
if the curvature of nonspherical particles at their contact points are included in a second- 
order an d higher-order analy sis, then hypostatic packings of such particles can indeed be 
jammed (IDonev et a/1 l2007al ). For example, ellipsoid packings are genera lly not jammed to 



2007al ) due to curvature 



first order in 5 but are jammed to second order in 5 (IDonev et al 
deviations from the sphere. 

To illustrate how nonspherical jammed packings can be hypostatic, Fig. [22] depicts two 
simple two-dimensional examples consisting of a few fixed ellipses and a central particle 
that is translationally and rotationally trapped by the fixed particles. Generically, four con- 
tacting particles are required to trap the central one. However, there are special correlated 
configurations that only require three contacting particles to trap the central one. In such 
instances, the normal vectors at the points of contact intersect at a common point, as is 
necessary to achieve torque balance. At first glance, such configurations might be dismissed 
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as probability-zero events. However, it was shown that such nongeneric configurations are 
degenerate (frequently encountered). This "focusing capacity" toward hypostatic values of 
Z applies to large jammed packings of nonspherical particles and in the case of ellipsoids 
must be present for s ufficiently small aspect ratios for a variety of realistic packing protocols 
( iDonev et all l2007al ). It has been suggested that the degree of nongenericity of the pack- 
ings be quantified by determining the fraction of local coordination configuratio ns in which 



the ce ntral particles have fewer contacting neighbors than the average value Z f lJiao et al. 
2010ah . 





FIG. 22 Simple examples of hypoconstrained packings in which all particles are fixed, except the 
central one. Left panel: Generically, four contacting particles are required to trap the central 
one. Right panel: Special correlated configurations only require three contacting particles to trap 
the central one. The normal vectors at the points of contact intersect at a common point, as is 
necessary to achieve torque balance. 



Having established that curvature deviatio ns from the spherica 



fundamental influence on constraint counting ( IDonev et al. 



reference shape exert a 



2007al ). it is clear that similar 



effects will emerge when the hard-particle interactions are replaced by nonspherical particles 
interacting with soft short-range repulsive potentials. It immediately follows that jamming 
to first and second order in S for hard nonspherical particles, for example, leads to quadratic 
and quartic modes in the vibrational energy spectrum for packings of such particles that 



inte ract with purely soft r epulsive interactions. The reader is referred to 



and 



Zeravcic et al 



Mailman et al. 



2009 



2009 



for studies of the latter type for ellipses and ellipsoids, respectively. 
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It is noteworthy that in striking contrast with MRJ-like sphere packings, the rattler 
concentrations of the MRJ-like ellipsoid p ackings appear pract ically to vanish outside of some 



small neighborhood of the sphere point (IDonev et al. 



2007c 



It was shown that MRJ-like 
packings of nearly-spherical ellipsoids can be obtained with <fi « 0.74, i.e., packing f raction s 
approaching those of the densest three-dimensional sphere packings (IDonev et all 12004a ). 
This suggested that there exist ordered ellipsoid packings with appreciably hi gher densities . 



Indeed , the densest known ellipsoid packings were subsequently discovered (IDonev et al. 



2004bl ); see Fig. [23j These represent a new family of non-Bravais lattice packings of ellipsoids 



with a packing fraction that always exceeds the density of the densest Bravais lattice packing 
(0 = 0.74048 . . .) with a maximal packing fraction of <fi = 0.7707. . ., for a wide range of 
aspect ratios (a < l/y/3 and a > \/3). In these densest known packings, each ellipsoid has 
14 contacting neighbors and there are two particles per fundamental cell. 

For identical ellipse packings in M 2 , the maximally dense arrangement is obtained by 
an affine stretching of the optimal tri angular-lattice packing of circula r disks with max = 



Til 



/vl2, which leaves 



unchanged (IDonev et al 



2004b 



Fejes Tothl . 11964 ): see Fig. E 



This maximally dense ellipse packing is not rotationally jammed for any noncircular shape, 



since it can 



) e shear ed continuously without introducing overlap or changing the density 



(IDonev et al. 



2007al ). The packing is, however, strictly translationally jammed. 



B. Superball Packings 

Virtually all systematic investigations of the densest particle packings have been carried 
out for convex objects. A c/-dimensional superball is a centrally symmetric body in d- 
dimensional Euclidean space occupying the region 

\xi\ 2p + \x 2 \ 2p + ■ ■ ■ + \x d \ 2p < 1, (29) 

where Xi (i = 1, . . . , d) are Cartesian coordinates and p > is the deformation parameter 
(not pressure as denoted in Sec. IIV.C1) . which controls the extent to which the particle shape 
has deformed from that of a rf-dimensional sphere (p = 1). Thus, superballs constitute a 
large family of both convex (p > 0.5) and concave (0 < p < 0.5) particles (see Fig. 126]) . 

In general, a "superdisk," the designation in the two-dimensional case, possesses square 
symmetry. As p moves away from unity, two families of superdisks with square symmetry 



55 



OA 06 0.6 1 L.2 1.4 1.6 IS 2 

Aspect ratio a 

FIG. 23 The packing fraction of the "laminated" non-Bravais lattice packing of ellipsoids (with a 
two-particle basis) as a function of the aspect ratio a. The point a = 1 corresponding to the face- 
centered cubic lattice sphere packing is shown, along with the two sharp maxima in the packing 
fraction for prolate ellipsoids with a = \/3 and oblate ellipsoids with a = l/y/3, as illustrated 
in the insets. For both a < l/y/3 and a > y/3, the packing fractions of the laminated packings 
drop off precipitously holding the particle orientations fixed (blue lines). The presently maximal 
achievable packing fraction eft = 0.7707 ... is highlighted with a thicker red line, and is constant for 
a < l/y/3 and a > v3 because there is an affi ne stretch by an arb itrary factor along a direction 



in a mirror plane of the particle directions; see 



Donev et al 



2004b 



can be obtained depending on whether p < 1 or p > 1 such that there is a 45-degree 
rotation with respect to the "protuberances." When p < 0.5, the superdisk is concave (see 
Fig. I25"j) . The candi date maximally d ense packings were recently proposed for all convex 



and concave shapes (jJiao et all 120081 ). These are achieved by two different families of 



Bravais lattice packings such that max is nonanalytic at the "circular-disk" point (p = 1) 
and increases significantly as p moves away from unity. The broken rotational symmetry 
of superdisks influences the packing characteristics in a non-trivial way that is distinctly 
different from ellipse packings. Recall that for ellipse packings, no improvement over the 
maximal circle packing density is possible. For superdisks, one can take advantage of the 
four-fold rotationally symmetric shape of the particle to obtain a substantial improvement 
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FIG. 24 A portion of the densest packing of congruent ellipses, which is simply an affine transfor- 
mation of the densest circle packing, i.e., the densest triangular lattice of circles. This maximally 
dense ellipse packing is not rotationally jammed for any noncircular shape, but is strictly transla- 
tionally jammed. 
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FIG. 25 (color online) Super disks with different values of the deformation parameter p. 



on the maximal circle packing density. By contrast, one needs to use higher- dimensional 
counterparts of ellipses (d > 3) in order to improve on max for spheres. Even for three- 
dimensional ellipsoids, (j)^^ incre ases smoothly as the aspect ratios of the semi-axes vary 
from unity ( iDonev et all l2004bl ). and hence has no cusp at the sphere point. In fact, 



congruent three-dimensional ellipsoid p ackings have a 



only when they are randomly jammed ( IDonev et all l2004al ) 



cusp- 



ike behavior at the sphere point 



Increasing the dimensionality of the particle imbues the optimal "superball" pa ckings with 



struc tural characteristics that are richer than their two-dimensional counterparts f lJiao et al. 



20091 ). For example, in three dimensions, a superball is a perfect sphere at p — 1, but can 
possess two types of shape anisotropy: cube-like shapes (three-dimensional analog of the 
square symmetry of the superdisk) and octahedron-like shapes, depending on the value of 
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FIG. 26 (color online). Superballs with different values of the deformation parameter p. 




(a) 



(b) 



(d) 



FIG. 27 (color online). Candidate optimal packings of superballs: (a) The Co-lattice packing of 
superballs with p = 1.8. (b) The Ci-lattice packing of superballs with p = 2.0. (c) The Oo-lattice 
packing of superballs with p = 0.8. (d) The ©i-lattice packing of superballs with p = 0.55. 



the deformation parameter p (see Fig.l26l). As p continuously increases from 1 to oo, we have 
a family of convex superballs with cube-like shapes; at the limit p = oo, the superball is a 
perfect cube. As p decreases from 1 to 0.5, a family of convex superballs with octahedron- 
like shapes are obtained; at p — 0.5, the superball becomes a regular octahedron. When 
p < 0.5, the superball still possesses an octahedron-like shape but is now concave, becoming 
a three-dimensional "cross" in the limit p — > 0. Note that the cube and regular octahedron 
(two of the five Platonic polyhedra) have the same group symmetry (i.e., they have the same 
48 space group elements) because they are dual to each other. Two polyhedra are dual to 
ea ch other if t he ve rtices of one correspond to the faces of the other. 

have obtained analytical constructions for the densest known super- 



Jiao et al 



2009 



ball packings for all convex and concave cases. The candidate maximally dense packings 
are certain families of Bravais lattice packings (in which each particle has 12 contacting 
neighbors) possessing the global symmetries that are consistent with the symmetries of a 
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super ball. Evidence is provided that these packings are indeed optimal, and Torquato and 



Jiao (ITorquato and Jiao . 



2009bl ) have conjectured that the densest packings of all convex 



superballs are their densest lattice packings; see Fig. [571 For superballs in the cubic regime 
(p > 1), the candidate optimal packings are achieved by two families of Bravais lattice 
packings (Co and Ci lattices) possessing two-fold and three-fold rotational symmetry, re- 
spectively, which can both be considered to be continuous deformations of the fee lattice. 
For superballs in the octahedral regime (0.5 < p < 1), there are also two families of Bra- 
vais lattices (Oq and Oi lattices) obtainable from continuous deformations of the fee lattice 



keeping its four-fold rotational symmetry, an c 



from the densest lattice packing for regular 



octahedra ( iBetke and Henkl . |2000| ; Minkowski 119051 1. keeping the translational symmetry of 
the projected lattice on the coordinate planes, which are apparently optimal in the vicinity 
of the sphere point and the octahedron point, respectively (see Fig. |27|) . 

The proposed maximal packing density max as a function of deformation parameter p is 
plotted in Fig. [28l As p increases from unity, the initial increase of max is linear in (p — 1) 
and subsequently max increases monotonically with p until it reaches unity as the particle 
shape becomes more like a cube, which is more efficient at filling space than a sphere. These 
characteristics stand in contrast to those of the densest known ellipsoid packings, achieved by 
certain crystal arrangements of congruent spheroids with a two-particle basis, whose packing 
density as a functio n of aspect ratios ha s zero initial slope and is bounded from above by a 



value of 0.7707 . . . (IDonev et al 



is linear in (1— p). Thus, 



2004bl ). As p decreases from unity, the initial increase of 
is a nonanalyti c function of p at p = 1, which is consistent 



2008 



). However, the behavior 



with conclusions made about superdisk packings (I Jiao et al . 
of 0max as the superball shape moves off the sphere point is distinctly different from that 
of optimal spheroid packings, for which </> max increases smoothly as t he aspect ratios of th e 
semi-axes vary from unity and hence has no cusp at the sphere point ( IDonev et all l2004bl ). 
The density of congruent ellipsoid packings (not <ft ma , y ) has a cusp-like beha vior at the sphere 



point only when the packings are randomly jammed (IDonev et all [2004aJ). The distinction 
between the two systems results from different broken rotational symmetries. For spheroids, 
the continuous rotational symmetry is only partially broken, i.e., spheroids still possess one 
rotationally symmetric axis; and the three coordinate directions are not equivalent, which 
facilitates dense non-Bravais packings. For superballs, the continuous rotational symmetry 
of a sphere is completely broken and the three coordinate directions are equivalently four-fold 
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FIG. 28 (color online) . Density versus deformation parameter p for the packings of convex super- 



balls (Jiao et al 



20091 ) . Insert: Around p* = 1.1509 . . ., the two curves are almost locally parallel 



to each other. 



rotationally symmetric directions of the particle. Thus, a superball is less symmetric but 
more isotropic than an ellipsoid, a shape characteristic which apparently favors dense Bravais 
lattice packings. The broken symmetry of superballs makes their shapes more efficient in 
tiling space and thus results in a larger and faster increase in the packing density as the 
shape moves away from the sphere point. 

As p decreases from 0.5, the superballs become concave particles, but they still possess 
octahedron-like shapes [see Fig. l2"6T a)]. The lack of simulation techniques to generate concave 
superball packings makes it very difficult to find the optimal packings for the entire range 
of concav e shapes (0 < P < 0.5). However, based on their conclusions for convex superball 



packings, 



Jiao et al 



2009 



conjectured that near the octahedron point, the optimal packings 
possess similar translational symmetry to that of the Oi-lattice packing, and based on 
theoretical considerations proposed candidate optimal packings for all concave cases with a 
de nsity versu s p as s hown in Fig{ 



Jiao et al 



2010al have determined the packing fractions of maximally random jammed 



(MRJ) packings of binary superdisks in M 2 and monodispersed superballs in R 3 . They found 
that the MRJ densities of such packings increase dramatically and nonanalytically as one 
moves away from the circular-disk and sphere point (p = 1). Moreover, these disordered 
packings were demonstrated to be hypostatic, i.e., the average number of contacting neigh- 
bors is less than twice the total number of degrees of freedom per particle, and the packings 
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FIG. 29 (c olor online) 



superballs (jjiao et al 



Dens ity versus deformation parameter p for the lattice packings of concave 



2009). Insert: a concave superball with p = 0.1, which will becomes a 



three-dimensional cross at the limit p — > 0. 



are mechanically stable. As a result, the local arrangements of the particles are necessarily 
nontrivially correlated to achieve jamming and hence "nongeneric." The degree of "non- 
genericity" of the packings was quantitatively characterized by determining the fraction of 
local coordination structures in which the central particles have fewer contacting neighbors 
than the average value Z. Figure [30] depicts local packing structures with more contacts 
than average and those with less contacts than average in MRJ binary superdisk packings 
for different p values. It was also explicitly shown that such seemingly "special" packing 
configurations are counterintuitively not rare. As the anisotropy of the particles increases, 
it was shown that the fraction of rattlers decreases while the minimal orientational order 
(as measured by the cubatic order metric) increases. These novel characteristics result from 
the unique rotational symmetry breaking manner of superdisk and superball particles. 



C. Polyhedron Packings 

Until some recent developments, very little was known about the densest packings of 
polyhedral particles. The difficulty in obtaining dense packings of polyhedra is related to 
their complex rotational degrees of freedom and to the non-smooth nature of their shapes. 

The Platonic solids (mentioned in Plato's Timaeus) are convex polyhedra with faces 
composed of congruent convex regular polygons. There are exactly five such solids: the 
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FIG. 30 (color online). Local packing structures with more contacts than average (shown in blue) 
and those with less contacts than average (shown in pink) in two-dimension al MRJ superdisk 



Jiao et al. 



2010a. 



packings for different values of the deformation parameter p, as obtained from 
Note that the degree of nongenericity decreases as p moves away from its sphere-point value. Here 
the size ratio (diameter of large superdisks divided by the diameter of the small superdisks) is 1.4 
and the molar ratio (number of large superdisks divided by the number of small superdisks) is 1/3. 



A 



P1 




FIG. 31 (color online). The five Platonic solids: tetrahedron (PI), icosahedron (P2), dodecahedron 
(P3), octahedron (P4), and cube (P5). 

tetrahedron (PI), icosahedron (P2), dodecahedron (P3), octahedron (P4), and cube (P5) 
(see Fig. I3T1) [We note in passi ng that viral capsids often have icosahedral symmetry; see, 



for example, 



Zandi et al. 



2004 ] Here we focus on the problem of the determination of the 
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densest packings of each of the Platonic solids in three-dimensional Euclidean space M 3 , 
except for the cube, which is the only Platonic solid that tiles space. 

It is useful to highlight some basic geometrical properties of the Platonic solids that we 
will employ in subsequent sections of this review. The dihedral angle 9 is the interior angle 
between any two face planes and is given by 

. 9 cos(7r/g) 



sin 



(30) 



2 sin(7r/p) ' 

where p is the number of sides of each face and q is the number of faces meeting at 



each vertex (ICoxeter 



19731 ). Thus, 9 is 2 sin \l/V3), 2sin- 1 ($/ v / 3), 2 sin" 1 ^/^ 



2 sin _1 (-\/2/3), and 7r/2, for the tetrahedron, icosahedron, dodecahedron, octahedron, and 
cube, respectively, where $ = (1 + \/5)/2 is the golden ratio. Since the dihedral angle for 
the cube is the only one that is a submultiple of 2ir, the cube is the only Platonic solid that 
tiles space. It is noteworthy that in addition to the regular tessellation of K 3 by cubes in the 
simple cubic lattice arrangement, there is an infinite number of other irregular tessellations 
of space by cubes. Th i s tiling -degeneracy example vividly illustrates a fundamental point 



made by 



Kansal et al. 



2002U namely, packing arrangements of nonoverlapping objects at 



some fixed density can exhibit a large variation in their degree of structural order. We note 
in passing that there are two regular dodecahedra that independently tile three-dimensional 
(negative l y cur ved) hyperbolic space H 3 , as well as one cube and one regular icosahedron 
( jCoxeterl . Il973l ); see Sec. Kill for additional remarks about packings in curved spaces. 

Every polyhedron has a dual polyhedron with faces and vertices interchanged. The dual 
of each Platonic solid is another Platonic solid, and therefore they can be arranged into dual 
pairs: the tetrahedron is self-dual (i.e., its dual is another tetrahedron), the icosahedron and 
dodecahedron form a dual pair, and the octahedron and cube form a dual pair. 

An Archimedean solid is a highly symmetric, semi-regular convex polyhedron composed 
of two or more types of regular polygons meeting in identical vertices. There are thirteen 
Archimedean solids: truncated tetrahedron (Al), truncated icosahedron (A2), snub cube 
(A3), snub dodecahedron (A4), rhombicosidodecahdron (A5), truncated icosidodecahdron 
(A6), truncated cuboctahedron (A7), icosidodecahedron (A8), rhombicuboctahedron (A9), 
truncated dodecahedron (A10), cuboctahedron (All), truncated cube (A12) and truncated 
octahedron (A13) (see Fig. I32|) . Note that the truncated octahedron is the only Archimedean 
solid that tiles R 3 . The duals of the Archimedean solids are an interesting set of new 
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A1 A2 A3 A4 A5 




A11 A12 A13 

FIG. 32 (color online) The 13 Archimedean solids: truncated tetrahedron (Al), truncated icosa- 
hedron (A2), snub cube (A3), snub dodecahedron (A4), rhombicosidodecahdron (A5), truncated 
icosidodecahedron (A6), truncated cuboctahedron (A7), icosidodecahedron (A8), rhomb icubocta- 
hedron (A9), truncated dodecahedron (A10), cuboctahedron (All), truncated cube (A12), and 
truncated octahedron (A13). This typical enumeration of the Archimedean solids does not count 
the chiral forms (not shown) of the snub cube (A3) and snub dodecahedron (A4), which implies 
that the left-handed and right-handed forms of each of these pairs lack central symmetry. The re- 
maining 11 Archimedean solids are non-chiral (i.e., each solid is superposable on its mirror image) 
and the only non-centrally symmetric one among these is the truncated tetrahedron. 

polyhedra that are called the Archimedean duals or the Catalan polyhedra, the faces of 
which are not regular polygons. 

Another important observation is that the tetrahedron (PI) and the truncated tetrahe- 
dron (Al) are the only Platonic and non-chiral Archimedean solids, respectively, that are 
not centrally symmetric. The chiral snub cube and chiral snub dodecahedron are the only 
other non-centrally symmetric Archimedean solids. A particle is centrally symmetric if it 
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has a center C that bisects every chord through C connecting any two boundary points of 
the particle, i.e., the center is a point of inversion symmetry. We will see that the central 
symmetry of the majority of the Platonic and Archimedean solids (P2 - P5, A2 - A13) dis- 
tinguish their dense packing arrangements from those of the non-centrally symmetric ones 
(PI and Al) in a fundamental way. 

Tetrahedral tilings of space underlie many different molecular systems 



( IConway and Torquatol . 120061 ). Since regular tetrahedra cannot tile space, it is of in- 



terest to determine the highest density that such packings of particles can achieve (one 
of Hilbert's 18th problem set). It is of interest to note that the densest Bravais-lattice 
packing of tetrahedra (which requires all of the tetrahedra to have the same orienta- 
tions) has <t> = 18/49 = .367... and each tetrahedron touches 14 others. Recently, 



Conway and Torquatol . 



2006 



showed that the maximally dense tetrahedron packing cannot 
be a Bravais lattice (because dense tetrahedron packings favor face to face contacts) and 
found non-Bravais lattice (periodic) packings of regular tetrahedra with « 0.72. One such 
packing is based upon the filling of "imaginary" icosahedra with the densest arrangement of 
20 tetrahedra and then arranging the imagi nary icosahedr a in th eir densest lattice packing 



configuration. Using "tetrahedral" dice, 



Chaikin et al 



20071 experimentally gene rated 



Jaoshvili et al. 



2010 



for a 



jammed disordered packings of such dice with rs 0.75; see also 
refined version of this work. However, because these dice are not perfect tetrahedra (vertices 
and edges are slightly rounded), a defi nitive conclu sion could not be reached. Using physical 



Chen 



2008 



discovered a dense periodic arrangement 



models and computer algebra system 
of tetrahedra with (ft = 0.7786 . . ., which exceeds the density of the densest sphere packing 
by an appreciable amount. 



Torquato and Jiao 



2009a b 



In an attempt to find even denser packings of tetrahedra 
have formulated the problem of generating dense packings of polyhedra within an adaptive 
fundamental cell subject to periodic boundary conditions as an optimization problem, which 
they call the Adaptive Shrinking Cell (ASC) scheme. Starting from a variety of initial 
unjammed configurations, this optimization procedure uses both a sequential search of the 
configurational space of the particles and the space of lattices via an adaptive fundamental 
cell that shrinks on average to obtain dense packings. This was used to obtain a tetrahedron 



packing consisting o 
( jTorquato and Jiao . 



72 pa rticles per fundamental cell with packing fraction = 0.782 



2009al ). Using 314 particles per fundamental cell and starting from an 
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"equilibrated" low-density liquid c onfiguration, the same aut 



packing fraction to 



0.823 . . . (ITorquato and Jiaol. 



rors were able to improve the 



interestingly lacks any long-range order. lHaji-Akbari et al 



2009b ). T his packing arrangement 



20091 numerically constructed 



a periodic packing of tetrahedra made of parallel stacks of "rings" around "pentagonal" 
di pyramids con s isting of 82 particles per fundamental cell and a density = 0.8503 . . .. 



Kallus et al. 



2010 



found a remarkably simple "uniform" packing of tetrahedra with high 



symmetry consisting of only four particles per fundamental cell (two "dimers" ) with packing 
fraction = ^ = 0.854700 .... A uniform packing has a symmetry (in this case a point 
inversion symmetry) that takes one tetrahedron to another. A dimer is composed of a 



2010a 



pair of regular tetrahedra that exactly share a common face. ITorquato and Jiao 
subsequently presented an analytical formulation to construct a three-parameter family of 
dense uniform dimer packings of tetrahedra again with four particles per fundamental cell. 
(A uniform dimer packing of tetrahedra has a point of inversion symmetry operation that 
takes any dimer into another.) Making an assumption about one of these parameters resulted 
in a two-parameter family, including those with a pack ing fraction as high as 4 



0.855506 . . . (see left panel of Fi g. El. 



was made in the formulation of 



Chen et al. 



12250 
14319 



Torquato and Jiao 



20 lOl re cognized that such an assumption 



2010al and employed a similar formalism 



to obtain a three-parameter family of tetrahedron packings, including the densest known 
dimer packings of tetrahedra with the very slightly higher packing fraction <fi = |^ = 



0.856347. . . (see right panel of Fig. 1331 . The most general analytical formulation to date 
to construct dense periodic packings of tetra hedra with four particles per fundamental cell 
was carried out by ITorquato and Jiaol . l2010bl . This study involved a six-parameter family of 
dense tetrahedron packings that includes as special cases all of the aforementioned "dimer" 
packings of tetrahedra, including the densest known packings with packing fraction <ft = 
= 0.856347 . . .. This most recent investigation strongly suggests that the latter set 
of packings are the densest among all packings with a four-particle basis. Whether these 
packings are the densest packings of tetrah edra among all packings is an open question for 



reasons given by 



Torquato and Jiaol . 



20101) 



Using the ASC scheme and a varie t y of in itial conditions with multiple particles in the 
fundamental cell, (ITorquato and Jiaol . |2009aJJb|) were also able to find the densest known 
packings of the octahedra, dodecahedra and icosahedra (three non-tiling Platonic solids) 
with densities 0.947 . . ., 0.904 . . ., and 0.836 . . ., respectively. Unlike the densest tetrahedron 
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FIG. 33 (color online) Closely related dense tetrahedron packings. Left panel: A portion of one member of 
the densest two-parameter family of tetrahedron packings wit h 4 particles per fundamental cell and packing 



fraction 



14319 



0.855506 . . . (ITorquato and Jiao 



2010al ). The blue particles (lighter shade) represent 



two dimers (i.e., 4 tetrahcdra) but from this perspective each dimer appears to be a single tetrahedron. 
Right panel: A portion of one member of the densest three-parameter family of tet rahedron packings 



with 4 particles per fund amental cell and packing fraction </> = = 0.856347... (jChen et al 



Torquato and Jiao 



2010; 



2010bf ). The green particles (lighter shade) represent two dimers. It is clear that these 



two packings are configurationally almost identical to one another. 



packing, which must be a non-Bravais lattice packing, the densest packings of the other non- 



tiling Platonic so 



lattice packings (IBetke and Henk 



ids found by the algor i thm are their pre viously known optimal (Bravais 



2000 



Minkowski 



19051 ): see Fig. 



These simulation 



results as well as other theoretical considerations, which we briefly describe immediately 
below, led them to general organizing principles concerning the densest packings of a class 
of nonspherical particles. 

Rigorous upper bounds on the maximal packing fraction max of packings of nonspherical 
particles of general shape can be used to assess the packing efficiency of a particular dense 
packing of such particles. However, it has been highly challenging to formulate upper bounds 
for non-tiling particle packings that are nontrivially less than unity. It has recently been 



shown that 



of a packing of congruent nonspherical particles of volume vp in IR 3 is 
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FIG. 34 (color online) Portions of the dens est lattice packings of thre e of the centrally symmetric 



2009 



. Left panel: Icosahedron 



Platonic solids found by the ASC scheme (jTorquato and Jiaol . 
packing with packing fraction eft = 0.8363. . .. Middle panel: Dodecahedron packing with packing 
fraction cj> = 0.9045 .... Right panel: Octahedron packing with packing fraction = 0.9473 .... 



bounded from above according to 



'max — (Pm&x 



mm 



Vp 71 



(31) 



where v$ is the volume of the largest sphere that can be inscribed in the nonspheri cal 



2009al Jbh. 



particle and 7r/Vl8 is the maximal sphere-packing density (jTorquato and Jiao 
The upper bound ( 1311) will be relatively tight for packings of nonspherical particles provided 
that the asphericity 7 (equal to the ratio of the circumradius to the inradius) of the particle 
is not large. Since bound ( 13TT) cannot generally be sharp (i.e., exact) for a non-tiling, 
nonspherical particle, any packing whose density is close to the upper bound (j31~j) is nearly 
optimal, if not optimal. It is noteworthy that a majority of the centrally symmetric Platonic 
and Archimedean solids have relatively small asphericities and explain the corresponding 
small differences between <>['.,.„. and the pa cking fraction of the densest lattice packing max 



flBetke and Henkl. 



2000 



Torquato and Jia 



iacl . 



Minkowski 



2009a 



19051 ). 



JbJ have demonstrated that substantial face-to-face contacts be- 
tween any of the centrally symmetric Platonic and Archimedean solids allow for a higher 
packing fraction. They also showed that central symmetry enables maximal face-to-face 
contacts when particles are aligned, which is consistent with the densest packing being the 
optimal lattice packing. 

The aforementioned simulation results, upper bound, and theoretical considerations led 
to the following three conjec tures concerning the dense st packin gs of polyhedra and other 



nonspherical particles in R 3 (jTorquato and Jiao 



200Mb, 



a, 



l2010bh : 



Conjecture 1: The densest packings of the centrally symmetric Platonic and Archimedean 
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solids are given by their corresponding optimal lattice packings. 

Conjecture 2: The densest packing of any convex, congruent polyhedron without central 
symmetry generally is not a (Bravais) lattice packing, i.e., set of such polyhedra whose 
optimal packing is not a lattice is overwhelmingly larger than the set whose optimal packing 
is a lattice. 

Conjecture 3: The densest packings of congruent, centrally symmetric particles that do 
not possesses three equivalent principle axes (e.g., ellipsoids) generally cannot be Bravais 
lattices. 

Conjecture 1 is the analog of Kepler's sphere conjecture for the centrally symmetric 
Platonic and Archimedean solids. Note that the densest known packing of the non-centrally 
symmetric truncated tetrahedron is a non-lattice packing with density at least as high as 



23/24 = 0.958333 . . . flConwav and Torquato . 



20061 ). The arguments leading to Conjecture 1 



also strongly suggest that the densest packing s of sup erballs are given by their cor responding 



optim al lattice packings (jTorquato and Jiaol . l2009bl ). which were proposed by 



Jiao et al. 



2009 



D. Additional Remarks 



: the Platonic and Archimedean 



It is noteworthy that the densest known packings of all o 
solids as well as the densest known packings of superballs (IJiao et all l2009f ) and ellipsoids 
( jDonev et all l2004bl ) in M 3 have packing fractions that exceed the optimal sphere packing 
value (p^ ax = tt/v^ = 0.7408 . . .. These results are consistent with a conjecture of Ulam 
who proposed without any justification [in a private communication to Martin Gardner 
( IGardneit l200lh ] that the optimal packing fraction for congruent sphere packings is smaller 
than that for any other convex body. The sphere is perfectly isotropic with an asphericity 7 
of unity, and therefore its rotational degrees of freedom are irrelevant in affecting its packing 
characteristics. On the other hand, each of the aforementioned convex nonspherical particles 
break the continuous rotational symmetry of the sphere and thus its broken symmetry can be 
exploited to yield the densest possible packings, which might be expected to exceed <j) max = 



7r/vT8 = 0.7408 . . . ( Torquato and Jiaol . |2009aj). However, broken rotational symmetry in 



and of itself may not be sufficient to satis 
little or no symmetry ( jTorquato and Jiao . 



y Ulam 's conjecture if the convex particle has a 



2009af ). 
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Apparently, the two-dimensional analog of Ulam's conjecture (optimal density of congru- 
ent circle packings (</> max = n/y/12 = 0.906899 . . .) is smaller than that for any other c onvex 
two-dimensional body) is false. The "smoothed" octagon constructed by (IReinhardtl . Il934l ) 
is conjectured to have smallest optimal packing fraction (0 max = (8— 4\/2— In 2) / (2\/2— 1) = 
0.902414 . . .) among all congruent centrally-symmetric planar particles. 

It will also be interesting to determine whether Conjecture 1 can be extended to other 
polyhedral packings. The infinite families of prisms and antiprisms provide such a class 
of packings. A prism is a polyhedron having bases that are parallel, congruent polygons 
and sides that are parallelograms. An antiprism is a polyhedron having bases that are 
parallel, congruent polygons and sides that are alternating bands of triangles. They are 
generally much less symmetric than either the Platonic or Archimedean solids. Moreover, 
even the centrally symmetric prisms and antiprisms generally do not possess three equivalent 
directions. Thus, it is less obvious whether Bravais lattices would still provide the optimal 



packing s for these solids, except for prisms that tile spac e , e.g., 



prisms (ITorquato and Jiaol . l2009bl ) . iTorquato and Jiao 



20091) 



rexagonal prism or rhombical 
have also commented on the 



validity of Conjecture 1 to polytopes in four and higher dimensions. 



XI. PACKING SPHERES IN HIGH-DIMENSIONAL EUCLIDEAN SPACES 



There has been resurgen t interest in spher e packings fo r d > 3 in both the p h ysica 



and mathematical sciences (lAdda-Bedia et al 



Cohn and Kumar 



2009 



Conway and Sloane 



2008 



1995 



Cohn 



Elkies 



2002 



2000 



Cohn and Elkies 



Frisch and Percusl. 



2003 



1999; 



Lue et al. 2010; van Meel et al. 2009; Parisi and Slanina. 2000; Parisi and Zamponi. 2006 



2010 



2006 



Rohrmann and Santos. 



2007 



Torquato and Stillinger 



Scardicchio et al. 



2006a 



2008 



van Meel et al. 



Skoge et al 



2006 



Torquato 



20091 1 . Remarkably, the optimal 



way of sending digital signals ov er noisy channels co r responds to the dense st sphere pack 



1998 



Shannon 



19481 ). These "error- 



ing in a high- dimensional space (jConway and Sloane . 
correcting" codes underlie a variety of systems in digital communications and storage, in- 
cluding compact disks, cell phones and the Internet; see Fig. [35j Physicists have studied 
sphere packings in high dimensions to gain insi ght into liquid and glassy states of matter as 



well a s phase behavio r in lower dimensions ([Adda-Bedia et al. 



well 


as 


1999 





van Meel et al. 



2009 



Parisi and Slanina 



2008 



Frisch and Percus 



2000l ; iParisi and Zamponi 



2006 



2010; 
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Rohrmann and Santosl . 120071 ; ISkoge et all 120061 ; Ivan Meel et all 120091 ) . Finding the densest 
packings in arbitrary dimen sion is a problem of long-standing interest in discrete geometry 



rtConway and Sloane 



19981 ) . A comprehensive review of this huge subject is beyond the 



scope of this article. We instead briefly summarize the relevant literature leading to a recent 



development that supports the counterint uitive possibility t 



for suf ficiently large d may be disordered (IScardicchio et al. 



r at th e densest sphere packings 



2008 



Torquato and Stillinger 



2006bl ). or at least possess fundamental cells whose size and structural complexity increase 
with d. 

The sphere packing problem seeks to answer the following question: Among all packings 
of congruent spheres in R d , what is the maximal packing density <t>-m^ and what are the cor- 



responding arrangements of the spheres (jConway and Sloane 



are known only for the first three space dimen sions ( 



known packings are Bravais lattice packings (jConway and Sloane 



Hale 



2005h. 



998l ) ? The optimal solutions 



or 4 < d < 9, the densest 



19981 ). For example, the 



"checkerboard" lattice D^, which is a d-dimensional generalization of the fee lattice (densest 
packing in R 3 ), is believed to be optimal in R 4 and R 5 . The remarkably symmetric E$ 
and Leech lattices in IR 8 and R 2 4 , res pectively, are most likely the densest packings in these 



20091 ). Table HTT1 lists the densest known sphere packings in 



dimensions flCohn and Kumar . 
R d for selected d. Interestingly, the non-lattice (periodic) packing P\q c (with 40 spheres 
per fundamental cell) is the densest known packing in R 10 , which is the lowest dimension 
in which the best known packing is not a (Bravais) lattice. It is noteworthy that for suffi- 
ciently large d, lattice packings are most likely not the densest (see Fig. [36]), but it becomes 
increasingly difficult to find explicit dense packing constructions as d increases. Indeed, the 
problem if finding the shortest lattice vector in a particular lattice packing (densest lattice 



packing) grows super-exponential 



polynomial-time hard) problems (lAjtai 



ly with d an d is in the class of NP-hard (non-deterministic 



19981 ). 



For la rge d, the best that one ca n do theoretically is to devise upper and lower bounds 



Minkowski 



1905 



on 4> max ( jConway and Sloand . 1199a ) . The nonconstructive lower bound of 
established the existence of reasonably dense lattice packings. He found that the maximal 
packing fraction </> max among all lattice packings for d > 2 satisfies 

((d) 



'max — 2<2-l 



(32) 



where ((d) = YlkLi k d is the Riemann zeta function. Note that for large values of d, the 
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FIG. 35 (color online) A fundamental problem in communications theory is to find the best way to 
send signals or "code words" (amplitudes at d different frequencies) over a noisy channel. Each code 
word s corresponds to a coordinate in M. d . The sender and receiver desire to design a "code book" 
that contains a large number of code words that can be transmitted with maximum reliability 
given the noise inherent in any communications channel. Due to the noise in the channel, the code 
word s will be received as r ^ s such that |s — r| < a, where the distance a is the maximum error 



associated with the channel. 



Shannon , 



19481 showed that the best way to send code words over 
a noisy channel is to design a code book that corresponds to the densest arrangement of spheres 
of radius 2a in M. d . To understand this remarkable result, this figure depicts two choices for the 
code book in two dimensions. Left panel: Three code words (depicted as green points) that are 
circumscribed by circles of radius a are shown in which the distance between any two words is 
chosen to be less than 2a. The received signal (red point) falls in the region common to all circles 
and therefore the receiver cannot determine which of the three code words shown was sent. Right 
panel: By choosing a code book in which the code words are at least a distance 2a apart (which 
corresponds to a sphere packing), any ambiguity about the transmitted code word is eliminated. 
Since one desires to send many code words per unit volume, the best code book corresponds to the 
densest sphere packing in R d in the limit that the number of code words tends to infinity. 



asymptotic behavior of the Minkowski lower bound is controlled by 2 



Since 1905, many exten s ions and generalization s of 



Conway and Sloane 



1998 



Davenport and Rogers! . 



1947 



have been o btained ( IB all 



1992 



Vance . 



20091 ). but none of these 
investigations have been able to improve upon the dominant exponential term 2~ d . It is 



useful to note that the packing fraction of a saturated packing of congruent spheres in 
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cT-1 




FIG. 36 (color online) Lattice packings in sufficiently high dimensio ns are almost surely unsatu rated 



1998 ). To 



because the "holes" or space exterior to the spheres dominates M rf (jConwav and Sloand . 
get an intuitive feeling for this phenomenon, it is instructive to examine the hypercubic lattice 
Z d (square lattice for M. 2 and the simple cubic lattice for IR 3 ). Left panel: A fundamental cell 
of Z d represented in two dimensions. Note that the distance between the point of intersection of 
the longest diagonal in the hypercube with the hypersphere boundary and the vertex of the cube 
along this diagonal is given by yd — 1 for a sphere of unit radius. This means that Z d already 
becomes unsaturated at d = 4. Placing an additional sphere in Z 4 doubles the density of Z 4 and, 
in particular, yields the four-dimensional checkerboard lattice packing D^, which is believed to 
be the optimal packing in R 4 . Right panel: A schematic "effective" distorted representation of 
the hypersphere within the hypercubic fundamental cell for large d, illustrating that the volume 
content of the hypersphere relative to the hypercube rapidly diminishes asymptotically. Indeed, 
the packing fraction of Z d is given by (f> = 7r rf / 2 /(r(l + d/2)2 d ). It is the presence of the gamma 
function, which grows like (d/2)l, in the denominator that makes Z rf far from optimal (except for 
d = 1). In fact, the checkerboard lattice with packing fraction (f> = 7r d / 2 /(r(l + d/2)2 ( - d+2 ^ 2 ) 
becomes suboptimal in relatively low dimensions because it too becomes dominated by larger and 
larger holes as d increases. 



for all d satisfies 



> 



2 d ' 



(33) 



The proof is trivial. A saturated packing of congruent spheres of unit diameter and packing 
fraction in M. d has the property that each point in space lies within a unit distance from 
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TABLE III The densest known sphere packings in IR rf for selected d. Except for the non-lattice 
packing Piq c in M 10 , all of the other densest known packings listed in this table are lattice packings: 
Z is the integer lattice, A2 is the triangular lattice, is the checkerboard lattice (a generalization 
of the fee lattice), Ed is one t he root lattices, and is the laminated lattice. The reader is referred 



to 



Conway and Sloand . 



1998 for further details. 



Dimension, d 


Packing Structure 


Packing Fraction, 


1 


Z 


1 


2 


A 2 


7r/Vl2 = 0.9068... 


3 


D 3 


n/y/l8 = 0.7404... 


4 


D 4 


tt 2 /16 = 0.6168. . . 


5 


D 5 


2tt 2 /(30 v / 2) = 0.4652... 


6 


E e 


3tt 2 /(144V3) = 0.3729... 


7 


E 7 


7r 3 /105 = 0.2952. . . 


8 


E$ 


vr 4 /384 = 0.2536 . . . 


9 


A 9 


27r 4 /(945V2) = 0.1457... 


10 


PlOc 


vr 5 /3072 = 0.09961... 


16 


Al6 


tt 8 /645120 = 0.01470. . . 


24 


A 24 


vr 12 /479001600 = 0.001929... 



the center of some sphere. Thus, a covering of the space is achieved if each sphere center is 
encompassed by a sphere of unit radius and the packing fraction of this covering is 2 d <ft > 1 . 
Thus, the bound ( 133|) . which is sometimes called the "greedy" lower bound, has the same 
dominant exponential term as ( 13~2"|) . 

We know that there exists a disordered but unsaturate d packing construction, known a s 
the "ghost" random sequential addition (RSA) packing (jTorquato and Stillingerl . l2006al ). 
that achieves the packing 2~ d for any d. This packing, depicted in Fig. [37] and described in 
its caption, is a generalization of the standard RSA packing, also described in the caption of 
Fig [37J It was shown that all of the n-particle correlation functions of this nonequilibrium 
model, in a certain limit, can be obtained analytically for all allowable densities and in 
any dimension. This represents the first exactly s olvable disordered sphere-packing model 



Matern 



19861 gave an expression for the pair correla- 



in arbitrary dimension. (Note that 
tion function for this model.) The existence of this unjammed disordered packing strongly 
suggests that Bravais-lattice packings (which are almost surely unsaturated for sufficiently 
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large d) are far from optimal for large d. Further support for this conclusion is the fact 
that the maximal "saturation" packing fraction of the standard d isordered RSA packin g 
apparently scales as d ■ 2~ d or possibly d ■ ln(d) • 2~ d for large d (ITorquato et all 120061 ) . 
Spheres in both the ghost and standard RSA packings cannot form interparticle contacts, 
which a ppears to be a crucial attribute to obtain exponential improvement on Minkowski's 



bound (ITorquato and Stillinger 



2006bl ). as we discuss below. 
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FIG. 37 (Color online) The ghost RSA packing model in R d is a subset of the Poisson point process 
and a generalization of the standard RSA process. The latter is produced by randomly, irreversibly, and 
sequentiall y placing nonoverlapping objects into a large volume in M. d that at some initial time is empty 



of spheres ([Cooper 



1988 



Feder 



1980: 



Viot et al 



1993 



Widom 



1966). If an attempt to add a sphere at 



some time results in an overlap with an existing sphere in the packing, the attempt is rejected and further 
attempts are made until it can be added without overlapping existing spheres. As the process continues in 
time, it becomes more difficult to find available regions into which the spheres can be added, and eventually 
in the saturation or infinite-time limit , no further additions arc possible. In the ghost RSA process, a 
sphere that is rejected is called a "ghost" sphere. No additional spheres can be added whether they overlap 
an existing sphere or a ghost sphere. The packing fraction at time t for spheres of unit diameter is given 
by 4>(t) = [1 — exp(— vi (l)t)]/2 d , where vi(R) is the volume of sphere of radius R [cf. J2J)]. Thus, we see 
that as t — > +oo, <j> — > 2~ d . This figure shows a configuration of 468 particles of a ghost RSA packing 
in a fundamental cell i n M 2 at a packing fraction very near its maximal value of 0.25, as adapted from 



Scardicchio et al 



2008. Note that the packing is clearly unsaturated and there are no contacting particles. 
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The best currently 



obtained by 



Ball . 



mown lower bound on 0^ax f° r dimensions not divisible by four was 



1992 



He found that 



> 



2(d-l)C(d) 
2 d 



(34) 



Vance 



2009 



For dimensions divisible by four 
on the maximal packing fraction among lattice packings: 

6d 



has recently found the tightest lower bound 



> 



(35) 



2 d e(l - e~ d )' 

Table [IV] gives the dominant asymptotic behavior of several lower bounds on 0^ax f° r large 
d. Note that the best lower bounds on </>n iax improve on Minkowski's bound by a linear 
factor in d rather than providing exponential improvement. This suggests that the packing 
fraction of the densest lattice packing in high d is controlled by the exponential factor 2~ d . 



TABLE IV Dominan t asymptotic behavior of lower bounds on 



d. The bound due to 



Vance 



for sphere packings for large 



20091 is applicable for dimensions divisible by four. The best lower 
bounds do not provide exponential improvement on Minkowski's lower bound; they instead only 
improve on the latter by a factor linear in d. 



(2)2- d 


Minkowski (1905) 


[ln( v / 2)d]2^ a! 


Davenport and Rogers (1947) 


(2d)2- d 


Ball (1992) 


(6d/e)2- d 


Vance (2009) 



TABLE V Dominant asymptotic behavior of upper bounds on 
d. 



for sphere packings for large 



(d/2)2-°- 5d 


Blichfeldt (1929) 


(d/e)2-°- 5d 


Rogers (1958) 


2-0.599CW 


Kabatiansky and Levenshtein (1979) 



Nontrivial upper bo unds on the max imal packing fraction max for any sphere packing in 



have been derived. 



Blichfeldt 



1929 



showed that the maximal packing fraction for all d 
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satisfies </> max < (of/2 + l)2" d / 2 . This upper bound was improved by iRogersl . Il958l . Il964l by 
an analysis of the Voronoi cells. For large d, R ogers' upper bound asymptotically becomes 



d-2^ 2 /e. 



Kabatiansky and Levenshtein 



1978 



found an even stronger bound, which in the 



limit d — > oo yields </> max < 2 _0 599M . All of these upper bounds prove that the maximal 
packing fraction tends to zero in the limit d — > oo. This rather counterintuitive high- 
dimensional property of sphere packings can be understood by recognizing that almost all of 
the volume of a <i-dimensional sphere for large d is concentrated near the sphere surface. For 
example, the volume contained with such a sphere up to 99% of it radius is (99/100) d , which 
tends to zero exponentially fast. Thus, in high- dimensional sphere packings (densest or not), 
almost all of the volume is occupied by the void space (space exterior to the spheres), which 
is to be contrasted with the densest sphere packings in low dimensions in which the volume 
co ntained within th e sphe res dominates over the void volume. 



Cohn and Elkies 



2003 



obtained and computed linear programming upper bounds, which 



provided improvement over Rogers' upper bounds for dimensions 4 through 36, but it is 



not yet 
large d. 



mown whether they imp rove upon the Kabatiaksky-Levenshtein upper bounds for 



Cohn and Kumar 



2009 



used these techniques to prove that the Leech lattice is the 
unique densest lattice in M 24 . They also proved that no sphere packing in IR 24 can exceed 
the density of the Leech lattice by a factor of more than 1 + 1.65 x 10~ 30 , and gave a new 
proof that E% is the unique densest lattice in R 8 . Table 53 provides the dominant asymptotic 
behavior of several upper bounds on </> max for large d. Note that the best upper and lower 
bounds on max differ by an exponential factor as d — > oo. 

Since 1905, many extens ions and generalizations of Minkowski's bound have been derived 



( iConway and S 



tial term 2 d . 



oane 



1998 0, but none o f them have improved upon the dominant exponen- 



Torquato and Stillingerl . l2006bl used a conjecture concerning the existence 



of disordered sphere packings and a ^-invariant optimization procedure that maximizes 
associated with a radial "test" pair correlation function <?2( r ) to provide the putative expo- 
nential i mprovement on Minkowski's 10 0-year-old bound on max . Specifically, a g2-invariant 
process (ITorquato and Stillingerl . l2Q02f ) is one in which the functional form of a "test" pair 
correlation ^(r) function remains invariant as density varies, for all r, over the range of 
packing fractions 

< <p < </>*. (36) 
The terminal packing fraction is the maximum achievable density for the ^-invariant 
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process subject to satisfaction of certain nonnegativity conditions on pair correlations. For 
any test g2(r) that is a function of radial distance r = |r| associated with a packing, i.e., 
<?2(j) — for r < D, they maximized the corresponding packing fraction, 



lim i 



subject to satisfying the following two necessary conditions: 



#2(7") > for all r, 



(37) 



(38) 



and 



S(k) > for all k. (39) 

Condition (|39|) is a necessary condition for the existence of any point process [cf. (J7J)]. When 
there exist sphere packings with a g 2 satisfying these conditions in the interval [0,0*], then 
one has the lower bound on the maximal packing fraction given by 



> 



(40) 



Torquato and Stillingerl . l2006bl conjectured that a test function (^(t") is a pair correlation 
function of a translationally invariant disordered sphere packing in M. d for < <ft < 0* for suffi- 
ciently large d if and only if the conditions ( )38l) and ( 1391) are satisfied. There is mounting evi- 
dence to support this conjecture. First, they identified a decorrelation principle, which states 
that unconstrained correlations in disordered sphere packings vanish asymptotically in high 
dimensions and that the g n for any n > 3 can be inferred entirely (up to small errors) from a 



knowledge of p and #2- T his decorrelation princip 
tion for the conjecture of iTorquato and Stillingerl. 



e, amo ng other results, provides justifica- 
, and is vividly ex hibited by the ex- 



2006b 



actly solvable ghost RS A packing process (ITorquato and Stillingerl . |2006aJ) as w ell as by com 



puter simulations in high dimensions of the max i mally random jammed state (ISkoge et al. 



on g 2 (jCostin and Lebowitz 



2004 



20061 ) and the standard RSA packing (iTorquatol. l2006r) . Second, other necessary condit ions 



Hopkins et al. 



20091 : 



Torquato and Stillinger 



2006b|) ap 



pear to only have relevance in very low dimensions. Third, one can recover the form of 
known rigorous bounds [cf. ( 13"2"|) and (I3~4"|) ] for specific test (72's when the conjecture is in- 
voked. Finally, in these two instances, configurations of disordered sphere packings on the 
torus have been numerically c onstructed with such go in low dimension s for densities up to 



the terminal packing fraction f lCrawford et al. 



2003 



Uche et al. 



20061 ) 



78 



Using a particular t e st pair correlation corresponding to a disordered sphere packing, 



Torquato and Stillinger 



2006bl found a conjectural lower bound on 



that is controlled 



by 2~(°- 77865 ---) d and the associated lower bound on the average contact (kissing) number 
Z is controlled by 2( 0,22134 ---)' i (a highly overconstrained situation). These results counter- 
intuitively suggest that the densest packings as d increases without bound may exhibit 
increasingly complex fundamental cells, or even become disordered at some sufficiently large 
d rather than periodic. The latter possibility would imply the existence of d isordered clas- 



Scardicchio et al. 



20081 demonstrated 



sical ground states for some continuous potentials, 
that there is a wide class of test functions (corresponding to disordered packings) that lead 
to precisely the same putative exponential improvement on Minkowski's lower bound and 
therefore the asymptotic form 2 _ (°- 77865 - ) q! is much more general and robust than previously 
surmised. 

Interestingly, the optimization problem define d abo ve is the dual of the infinite- 



dimensional linear program ( LP) devised by 



Cohnl. 



maximal packing fraction; see ICohn and Elkies 



2003 



2002 



to obtain upper bounds on the 



for a proof. In particular, let f(r) be 



a radial function in IR d such that 



/(r) < for r > D, 
f(k) > for all k, 



(41) 



where f(k) is the Fourier transform of f{r). Then the number density p is bounded from 
above by 

minJ^L. (42) 
2 d /(0) 

The radial function f(r) can be physically interpreted to be a pair potential. The fact 



that its Fourier transform must be nonnegative for all k is a wel 



dition for many-particle systems with pairwise interactions (IRuelle 



mown stability con- 



1999|). We see that 



whereas the LP problem specified by (1381) and (1391) utilizes information about pair correla- 
tions, its dual program fill) and fi2|) uses information about pair interactions. As noted by 



Torquato and Stillinger 



2006bl . even if there does not exist a sphere packing with g2 satis- 



fying conditions ( |38i) . ( 1391) and the hard-core constraint on g2, the terminal packing fraction 
0* can never exceed the Cohn-Elkies upper bound. Every LP has a dual program and when 
an optimal solution exis ts, there is no dual i ty gap between the upper bound and lower bound 



formulations. Recently, 



Cohn and Kumar 



2007al proved that there is no duality gap. 
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XII. REMARKS ON PACKING PROBLEMS IN NON-EUCLIDEAN SPACES 



Particle packing problems in non- Euclidean ( curved) spaces have been the focus of re- 



searc h in a vari ety of fiel d s, inc l uding physics (IBowick et al 



20071). biology ( 



Torquato et al. 



Goldberg 



2002; 



1967 



200G 



Prusinkiewicz and Lindenmayer 



Zandi et al 



Modes and Kamien . 



199( 



2004 ). co mmunications theory (IConwav and Sloane . 



1998 ). and geometry (ICohn and Kumar . 



2007b 



Conway and Sloane . 



(]~ Tamme J 



1930 



1998 



Hardin and Safn . 



20041 ) . Although a comprehensive overview of this topic is beyond the scope of this review, 



we highlight here some of the developments for the interested reader in spaces with constant 
positive and negative curvatures. We will limit the discussion to packing spheres on the 
positively curved unit sphere C M d and in negatively curved hyperbolic space M d . 

The kissing (or contact) number r is the number of spheres of unit radius that can 
simultaneously touch a unit sphere S' d ~ 1 (j Conway and Sloand . Il998l ). The kissing number 



problem asks for the maximal kissing number r max in R d . The determination of the maximal 
kissing number in M 3 spurred a famous debate between Issac Newton and David Gregory 
in 1694. The former correctly thought the answer was 12, but the latter wrongly believed 
that 13 unit spheres could simultaneously contact another unit spher e. The optim al kissing 
num ber T max in dimensions greater than three is o nly known for R 4 (jMusinl . 120081 ) , M 8 and 



d21 



(ILevenshteinl . 



1979 



Odlyzko and Sloand . 



19791 ) . Table IVII lists the largest known kissing 



numbers in selected dimensions. 

In geometry and coding theory, a spherical code with parameters (d, N, t) is a set of 
iV points on the unit sphere S^ -1 such that no two distinct points in that set have inner 
product greater than or equal to t, i.e., the angles between them are all at least cos -1 1. The 
fundamental problem is to maximize N for a given value of t, or equivalently to minimize t 
given N [sometimes cal led the Tammes problem, which was motivated by an application in 



bo tany ( 



by 



Tammes . 



193ol)1 



One of 



; he first rigorous studies of s pherical codes was presented 



Schuette and van der Waerden 



1951 



Delsarte et al 



1977 introduced much of the most 



important mathematical machinery to understand spherical codes and designs. One natural 
generalization of the best way to distribute points on S - 1 (or R d ) is the energy minimization 
problem: given some potential function depending on the pairwise distances between points, 
how should the points be arranged so as to minimize the total energy (or what are the 
ground-state configurations)? The original Thomson problem of "spherical crystallography" 
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TABLE VI The largest known kissing numbers for identical spheres in M. d for selected d. Except 
for M. 9 and IR 10 , the largest known kissing numbers listed in this table are those found in the densest 
lattice packings listed in Table IIII1 Pg a is a non-lattice packing with an average kissing number of 
235 3/5 but with a maximum kissing number of 306. P±ob is a non-lattice packing with an average 
kissing number of 340 1/3 but with a maximum kissing number of 500. The reader is referred to 



Conway and Sloane 



1998 for further details. 



Dimension, d 


Packing Structure 


Kissing Number, r 


1 


Z 


2 


2 


A 2 


6 


3 


D 3 


12 


4 




24 


5 


D, 


40 


6 


E 6 


72 


7 


E 7 


126 


8 


E 8 


240 


9 


P9a 


306 


10 


PlOb 


500 


16 


Aie 


4,320 


24 


A 24 


196,560 



seeks the ground states of electron shells interacting via the Coulomb potential; but it 
is als o profitable to stud y ground states of particle s interacting with other potentials on 



S^ 1 flBowick et all 



2006). 



Cohn and Kumar 



2007b have introduced the beautiful idea of a 



universally optimal configuration, a unique configuration that minimizes a class of potentials. 
In particular, they proved that for any fixed number of points N on S d ~ l there is a universally 
optimal configuration that minimizes all completely monotonic potential functions (e.g., all 
inverse power laws). 

Th e optimal spherical cod e problem is related to the densest local packing (DLP) problem 
in R d (jHopkins et all . l2010al ). which involves the placement of N nonover lapping spheres of 
unit diameter near an additional fixed unit-diameter sphere such that the greatest radius 
R from the center of the fixed sphere to the centers of any of the N surrounding spheres is 
minimized. Let us recast the optimal spherical code problem as the placement of the centers 
of N nonoverlapping spheres of unit diameter onto the surface of a sphere of radius R such 
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that R is minimized. It is has been proved that for any d, all solutions for R between unity 
and the golden ratio r = (1 + v5)/2 to the optimal spherical code proble m for N spheres are 
also solutions to the corresponding DLP problem (IHopkins et a/.l . l2010bl ). It follows that for 
any packing of nonoverlapping spheres of unit diameter, a spherical region of radius R less 
than or equal to r centered on an arbitrary sphere center cannot enclose a number of sphere 
centers greater than one more than the number that than can be placed on the spherical 
region's surface. 

We saw in Sec JVIHl that monodisperse circle (circular disk) packings in M 2 have a great 
tendency to crystallize at high densities due to a lack of geometrical frustration. The hyper- 
bolic plane H 2 (for a particular constant negative curvature, which measures the deviation 
from the flat Euclidean plane) provides a two-dimensional space in which global crystalline 
order in dense circle packings is frustrated, and thus affords a means to use circle pack- 
ings to understand fundame n tal fe atures of simple liquids, disordered jammed states and 



glasses. 



Modes and Kamien 



20071 formulated an expression for the equation of state for 
disordered hard disks in H 2 and compared it to corresp onding results obtained from molec- 



Modes and Kamien 



20071 derived a generalization of the virial 



ular dynamics simulations, 
equation in H 2 relating the pressure to the pair correlation function and developed the ap- 
propriate setting for extending integral-equation approaches of liquid-state theory. For a 
discussion of the mathematical subtleties associated with fin ding the densest packing s of 
identical ci-dimensional spheres in H d , the reader is referred to 



Bowen and Radin 
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XIII. CHALLENGES AND OPEN QUESTIONS 

The geometric-structure approach advanced and explored in this review provides a com- 
prehensive methodology to analyze and compare jammed disk and sphere packings across 
their infinitely rich variety. This approach also highlights aspects of present ignorance, thus 
generating many challenges and open questions for future investigation. Even for identical 
spheres, detailed characterization of jammed structures across the simple two-dimensional 
((f>-ip) order maps outlined in Sec. |V] is still very incomplete. A partial list of open and 
challenging questions in the case of sphere packings includes the following: 



1. Are the strictly jammed "tunneled" crystals ( jTorquato and Stillingerl . 120071 ) the family 



of lowest density collectively jammed packings under periodic boundary conditions? 
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2. How can the extremal jammed packings that inhabit the upper and lower boundaries 
of occupied regions of each of those order maps be unambiguously identified? 

3. What would be the shapes of analogous occupied regions if the two-parameter versions 
illustrated in Fig. [TBI were to be generalized to three or more parameters? 

4. To what extent can the rattler concentration in collectively or strictly sphere packings 
be treated as an independent variable? What is the upper limit to attainable rattler 
concentrations under periodic boundary conditions? 

5. What relations can be established between order metrics and geometry of the corre- 
sponding configurational-space polytopes? 

6. Can upper and lower bounds be established for the number of collectively and/or 
strictly jammed states for TV spheres? 

7. Upon extending the geometric-structure approach to Euclidean dimensions greater 
than three, do crystalline arrangements with arbitrarily large unit cells or even disor- 
dered jammed packings ever provide the highest attainable densities? 

Jamming characteristics of nonspherical and even non-convex hard particles is an area of 
research that is still largely undeveloped and therefore deserves intense research attention. 
Many of the same open questions identified above for sphere packings are equally relevant 
to packings of nonspherical particles. An incomplete list of open and challenging questions 
for such particle packings includes the following: 

1. What are the appropriate generalizations of the jamming categories for packings of 
nonspherical particles? 

2. Can one devise incisive order metrics for packings of nonspherical particles as well as 
a wide class of many-particle systems (e.g., molecular, biological, cosmological, and 
ecological structures)? 

3. Can sufficient progress be made to answer the first two questions that would lead to 
useful order maps? If so, how do the basic features of order maps depend on the shape 
of the particle? What are the lowest density jammed states? 
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4. Why does the deformation of spheres into ellipsoids have the effect of strongly inhibit- 
ing (if not completely preventing) the appearance of rattlers in disordered ellipsoid 
packings for a wide range of aspect ratios? Does this same inhibitory phenomenon 
extend to disordered jammed packings of other nonspherical shapes? 

5. For non-tiling nonspherical particles, can an upper bound on the maximal packing 
fraction be derived that is always strictly less than unity? (Note that upper bound 
(13 ip is strictly less than unity provided that the asphericity is sufficiently small.) 

Jamming characteristics of particle packings in non- Euclidean spaces is an area of research 
that is still largely undeveloped and therefore deserves research attention. This is true even 
for simple particle shapes, such as spheres. For example, we know very little about the 
jamming categories of spheres in curved spaces. Does curvature facilitate jamming or not? 
Does it induce ordering or disordering? These are just a few of the many challenging issues 
that merit further investigation. 

In view of the wide interest in packing problems that the research community displays 
there is reason to be optimistic that substantial conceptual advances are forthcoming. It 
will be fascinating to see how future review articles covering this subject document those 
advances. 
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